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SECTION  I 
INTRODUCTION 

This  report  covers  methods  based  on  an  optimality  criterion  to  design  a 
minimum  weight  structure.  These  methods  can  be  called  indirect  methods, 
since  the  objective  is  to  obtain  a  design  that  satisfies  a  certain  specified 
criterion  and  by  so  doing  indirectly  minimize  the  weight  of  the  structure. 

The  criterion  may  be  intuitive  or  derived  mathematically  based  on  the  nature 
of  the  problem.  We  will  concentrate  on  a  mathematical  criterion  derived 
by  differentiating  the  Lagrangian  with  respect  to  the  design  variables. 

In  the  optimality  criterion  methods  the  criterion  is  derived  for  the 
dominant  type  of  constraint  imposed  on  the  structure,  and  that  criterion 
is  used  to  develop  the  algorithm.  The  methods  are  iterative  in  nature 
because  of  the  nonlinearity  of  the  constraints  and  the  statical  indeterminacy 
of  the  structure.  In  deriving  the  optimality  criterion  and  developing  the 
algorithm,  full  use  is  made  of  the  knowledge  of  the  behavior  of  the  con¬ 
straints  imposed  on  the  structure.  The  algorithms  are  efficient  because 
they  are  specifically  developed  for  structural  optimization  and  generally 
treat  the  constraints  that  weakly  affect  the  behavior  of  the  structure  as 
passive  constraints. 

The  early  discussion  on  the  optimality  criterion  methods  as  applied 
to  a  discretized  structure  for  different  types  of  constraints  may  be  found 
in  References  1  through  4.  The “appl ication  of  the  method  to  the  stress- 
displacement-constrained  problem,  with  a  variation  in  the  recurrence 
relations  and  the  methods  used  to  estimate  the  Lagrange  multipliers  is 
discussed  in  References  5  through  20.  The  general  instability  of  the 
structure  as  a  constraint  on  the  structure  is  considered  in  References  21 
and  22.  The  optimization  of  a  structure  with  dynamic  loads  is  discussed 
in  References  23  through  25.  These  are  a  few  of  the  references  in  the 
general  field  of  structural  optimization  based  on  the  optimality  criterion 
approach.  Additional  references  may  be  found  in  the  recent  survey  paper 
(Reference  26),  and  a  detailed  review  of  the  state-of-the-art  in  Reference 
27.  In  Reference  28  a  comparison  and  the  relationship  between  the  dif¬ 
ferent  optimality  criterion  based  algorithms  are  presented.  The  relation¬ 
ship  between  the  optimality  criterion  and  mathematical  programming  methods 
is  presented  in  Reference  29.  The  relationship  between  the  structural 
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algorithm;,  ha  .vii  on  m  optimality  criterion  and  the  projection  method  is 
discussed  in  References  30  and  31.  Some  of  the  computer  programs  based 
on  the  optimality  criterion  to  design  minimum  weight  structures  are  docu¬ 
mented  in  References  32  through  36. 

The  constraints  imposed  on  the  structure  may  include  the  maximum 
allowable  stress  in  each  element,  the  displacement  limits  at  one  or  more 
locations,  system  stability,  dynamic  stiffness,  local  element  buckling 
etc.  In  addition  to  these  there  would  be  limitations  on  the  minimum  and 
maximum  sizes  of  the  elements.  "An  optimality  criterion  can  be  derived 
that  includes  all  these  constraints,  and  it  may  be  desirable  to  find  a 
design  that  satisfies  this  criterion."  However,  to  develop  an  efficient 
algorithm  based  on  such  a  criterion  and  effectively  handle  all  types  of 
constraints  would  be  impractical  and  generally  unnecessary.  For  most 
problems  it  is  developing  the  algorithm  that  is  more  difficult  than 
deriving  the  optimality  criterion.  In  the  case  of  most  structures  it  is 
likely  that  one  can  predict  the  type  of  constraint  which  will  be  the  most 
active  at  the  optimum  and  use  the  algorithm  based  on  that  constraint. 

Then  one  can  treat  all  other  constraints  as  passive  constraints.  It  is 
highly  unlikely  that  all  types  of  constraints  will  be  active  at  the 
optimum.  Sometimes  this  point  of  view  may  not  be  correct  and  will  not 
give  us  an  absolute  minimum  weight  design.  However,  it  will  give  us  a 
near  minimum  weight  design,  and  the  corresponding  optimization  algorithm 
will  be  efficient  and  easy  to  use  for  a  structure  with  a  large  number  of 
design  variables.  Even  with  this  approximation  to  the  overall  problem, 
when  the  total  number  of  constraints  of  the  same  type  are  large,  as  we 
will  see,  it  is  advantageous  to  make  additional  approximations  to  reduce 
the  computational  effort. 

The  optimality  criterion  derived  for  all  the  constraints  imposed  on 
the  structure  is  equivalent  to  the  Kuhn-Tucker  conditions  of  nonlinear 
mathematical  programming.  However,  in  deriving  the  optimality  criterion 
and  the  corresponding  structural  optimization  algorithm,  some  of  the  con¬ 
straints  are  treated  as  side  constraints  in  order  to  simplify  the  algorithm. 
A  good  example  of  this  is  the  minimum  and  maximum  size  limits  on  the  design 
variables.  These  constraints  generally  are  not  included  in  the  constraint 
equations  and  do  not  enter  the  optimality  criterion.  The  optimality 
criterion  derived  for  structural  optimization  for  a  particular  type  of 
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constraint  gives  information  on  the  distribution  of  energy  in  the 
structure  necessary  to  have  a  minimum  weight  design.  The  nature  of  the 
energy  depends  upon  the  type  of  constraint. 

The  optimization  procedure  may  be  divided  into  two  major  steps. 

These  are  (1)  the  analysis  of  the  structure,  and  (2)  the  redistribution  of 
the  material  so  that  the  weight  of  the  structure  is  reduced  and  the  active 
constraints  are  satisfied.  These  two  steps  are  repeated  alternately  until 
a  satisfactory  design  is  obtained.  The  analysis  of  the  discretized 
structure  is  performed  by  the  finite  element  method.  The  redistribution 
of  the  material  is  carried  out  by  using  a  recurrence  relation.  The 
recurrence  relation  modifies  the  design  variables  so  that  in  the  design 
space  the  initial  design  is  moved  towards  a  design  that  satisfies  the 
optimality  criterion.  The  recurrence  relation  contains  two  sets  of 
unknown  terms.  The  first  set  is  related  to  the  gradient  of  the  constraints 
and  the  second  set  is  the  Lagrange  multipliers.  It  is  necessary  to  deter¬ 
mine  these  unknowns  before  the  recurrence  relation  can  be  used. 

The  efficiency  and  the  convergence  behavior  of  the  algorithm  depend 
on  (1)  the  recurrence  relation  used  to  modify  the  design  variables, 

(2)  the  nature  of  the  approximations  made  to  derive  the  mathematical 
expression  for  the  unknowns  in  the  recurrence  relations,  and  (3)  how  these 
unknowns  are  determined. 

In  deriving  the  optimality  criterion  and  the  algorithm  based  on  it, 
we  will  assume  that  the  structure  is  discretized  into  a  number  of  elements, 
and  that  it  is  suitable  for  finite  element  analysis.  In  addition  we  will 
assume  that  the  only  design  variables  are  the  cross-sectional  areas  of 
the  bar  elements  and  the  thicknesses  of  the  plate  elements.  The  geometry 
of  the  structure,  materials,  the  loads  applied  to  it  are  fixed. 

This  report  discusses  different  methods  developed  to  optimize  struc¬ 
tures  with  constraints  on  displacements,  stresses,  minimum  sizes  and 
general  instability.  The  detailed  discussion  is  primarily  devi. ced  to 
constraints  on  stresses  an*  displacements,  since  these  are  the  basic 
constraints  in  any  structure.  First,  the  equations  of  the  displacement 
method  of  finite  element  analysis  are  reviewed;  then  the  optimization 
methods  are  presented.  Also  discussed  is  ti  rcect  of  the  reciprocal 
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design  variable  on  the  optimality  criterion  and  the  optimization  algorithms. 
The  general  optimization  procedure  and  the  different  approaches  that  one 
can  use  to  optimize  a  structure  are  then  presented.  The  report  concludes 
with  illustrative  problems  indicating  the  effect  of  the  different  algorithms, 
then  a  summary,  then  an  appendix.  The  appendix  contains  details  of  a 
three-bar  truss  problem  solved  by  using  different  algorithms. 
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SECTION  II 

BASIC  EQUATIONS  OF  ANALYSIS 

DISPLACEMENT  METHOD 

The  equations  of  finite  element  analysis  which  we  will  use  to  derive 
the  optimality  criterion  and  the  optimization  algorithm  are  summarized  in 
this  section. 

The  load  vector  { P>  is  related  to  the  displacement  vector  (u)  by 
the  equation 

[K]  {u>  =  {P}  (1) 

where  [K]  is  the  global  stiffness  matrix  of  the  structure.  The  strain 
energy  stored  in  the  ith  element  is  given  by 

Vi  =  (ul*  [k].  {u>.  (2) 

where  { u} -  is  the  displacement  vector  associated  with  the  ith  element  and 
[k].  is  the  ith  element  stiffness  matrix.  The  ith  element  strain  vector 
{e } -  is  related  to  {u}^  by 

te).  =  [B]i  {u}.  (3) 

where  [B]^  is  the  strain-displacement  matrix  associated  with  the  ith 
element. 

The  relation  between  the  stress  vector  {o>i  and  the  strain  vector 
{ e > -  in  the  ith  element  is  given  by 

<°}.  =  [E].  {el.  (4) 

where  [E]^  is  the  stress-strain  matrix  of  the  elastic  constants. 

Using  Eqs.  3  and  4,  a  relation  between  {o}^  and  { u> ^  can  be  written 
as 

io\  --  [D]i  (uli 

where  [D]^  is  the  strec<f  disp.acement  matrix  and  is  given  by 


(5) 
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by 


The  displacement  u.  at  the  jth  node  point  in  the  structure  is  given 

J 


Uj 


n  Q.  . 
v  _LL 
L  a 
i  =1  fti 


(7) 


or 


Uj 


Z  Q.  .Z. 
1-1  1 


(8) 


where  n  is  the  number  of  elements  and  zi  is  the  reciprocal  design  variable 
related  to  the  direct  design  variable  A.,  by 


Zi  =  A. 


(9) 


In  Eqs .  7  and  8,  Q^.  is  the  flexibility  coefficient  and  is  given  by 


(10) 


or 


If  j  "  1; 


(ID 


where  {sJ>.  is  the  virtual-displacement  vector  associated  with  the 

*  j 

virtual -load  vector  { S°> .  The  displacement  u.  is  also  given  by 

J 

u.  =  (u> t{Sj) 

J 


(12) 


For  a  bar  structure,  the  deflection  at  a  point  is  given  by  the  well-known 
relation 


_  "  TiVi 

J'  i=l  Vi 


(13) 


and  the  flexibility  coefficient  can  be  written  as 


Et 


(14) 
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where  T.  and  t?  are  the  forces  in  the  ith  bar  due  to  the  applied  load  { P } 
and  the  virtual  load  (SJ)  respectively. 

The  linear  stability  of  the  structure  is  determined  by  solving  the 
eigenvalue  problem 

[[K]  +  y[KG]]{n}  =  0  (15) 


where  [Kg]  is  the  geometric  stiffness  matrix  of  the  structure,  and  p  is 
the  eigenvalue.  The  eigenvalue  p.  is  given  by 

J 


(n}j[K]{n}j 

(n)  [KgKnlj 


(16) 


where  { n }  -  is  the  buckling  mode.  The  buckling  load  of  the  structure 

J  _ 

is  given  by  p{P},  where  p  is  the  lowest  eigenvalue. 

SCALING  OF  THE  DESIGN 

In  an  optimization  algorithm,  it  is  convenient  to  obtain  a  feasible 
design  after  each  iteration.  This  helps  keep  track  of  the  reduction  in 
the  weight  of  the  structure  after  each  iteration  and  also  helps  to  pick 
the  most  active  constraints.  A  feasible  design  can  be  obtained  by  scaling 
the  design  to  satisfy  the  specified  constraints.  The  design  variable  A^ 
can  be  written  as 

At  =  Ac.  (17) 

where  A  is  the  scaling  parameter.  In  Eq.  17,  a.  is  the  relative  value  of 
the  design  variable  A^ .  The  scaling  parameter  A  is  the  same  for  all  the 
elements.  Using  the  definition  of  A.,  given  in  Eq.  17,  Eq.  1  can  be 
written  as 


A[K' ]{ u)  =  { P>  (18) 

where  [ K ' ]  is  the  stiffness  matrix  of  the  structure  obtained  by  using  the 
relative  design  vector  (a).  If  { u ' >  is  the  relative  displacement  vector, 
then  the  equilibrium  equation  can  be  written  also  as 

[ K '  ] { u ' }  =  (P)  (19) 
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Comparing  Eqs.  18  and  19,  the  relationship  between  the  actual  displace¬ 
ment  vector  {u}  and  the  relative  displacement  vector  { u ' }  can  be  written 
as 

{u1 }  =  A  { u }  (20) 

Similarly,  expressing  the  strain-displacement  relation  and  the  stress- 
strain  relation  in  terms  of  the  relative  design  variable  {a},  the  relation¬ 
ship  between  the  relative  strains  and  stresses  can  be  obtained.  This  gives 


{e’}i  =  AU). 

(21) 

{a')i  =  A{a } i 

(22) 

where  { e 1 } ^  and  { a ' } ^  are  the  relative  strains  and  stresses  associated 
with  the  ith  element.  Eqs.  20  through  22  show  that  one  can  select  the 
parameter  A  to  satisfy  the  displacement  constraint  at  either  a  node  point 
or  the  stress  constraint  in  an  element.  The  A  corresponding  to  the  most 
critical  constraint  also  will  satisfy  all  other  constraints.  The  design 
with  the  actual  values  of  the  design  variables  satisfying  all  constraints 
is  given  by  Eq.  17,  where  A  is  determined  for  the  most  critical  constraint. 
The  simple  scaling  procedure  discussed  here  can  be  used  when  the  stiffness 
matrix  is  a  linear  function  of  the  design  variable.  In  other  cases,  one 
must  iterate  to  obtain  an  acceptable  scaling  parameter.  For  certain 
problems,  the  scaling  procedure  may  not  be  employable. 

If  the  scaling  parameter  A  is  Introduced  in  the  stability  equation 
(Eq.  15),  we  obtain 

[[K'J  +  u'CK^Kn’}  -  0  (23) 

where 


and 

(o')  =  A(n)  (25) 


Using  Eq.  24,  we  can  scale  the  design  so  that  the  buckling  load  of  the 
structure  is  equal  to  the  applied  load. 
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SECTION  III 

DISPLACEMENT  CONSTRAINTS 

The  displacement  constraints  are  fundamental  to  the  development  of 
an  optimization  algorithm  that  uses  the  displacement  method  of  finite 
element  analysis.  The  algorithms  for  other  types  of  constraints  are 
similar  to  that  of  the  displacement-constrained  problem.  In  the  case  of 
a  stress-constrained  problem,  for  example,  an  effective  algorithm  can  be 
developed  only  by  replacing  the  stress  constraints  by  the  equivalent  dis¬ 
placement  constraints.  We  will  first  consider  the  problem  of  designing 
a  structure  with  a  single-displacement  constraint;  then  discuss  the 
multiple  constraints.  The  index  1 i1  in  this  section  and  subsequent 
sections  refers  to  an  element  and  goes  from  1  to  n  where  n  is  the  number 
of  elements. 

SINGLE-DISPLACEMENT  CONSTRAINT 
Optimality  Criterion 

The  weight  of  the  structure  can  be  written  as 

n 

W(A)  =  E  p.-A.Jl,  (26) 

i=l  1  1  1 

where  is  the  specific  weight.  A.,  is  the  design  variable,  and  i.  is  a 
function  of  the  geometry  of  an  element.  The  product  A^  represents  the 
volume  of  the  element.  For  a  bar  element,  A-  is  the  cross-sectional  area 
and,  for  a  plate  element,  it  is  the  thickness.  Eq.  26  assumes  that  the 
weight  of  the  structure  is  a  linear  function  of  the  design  variable  Ai . 

A  single-displacement  constraint  can  be  written  as: 

91  =  C,  -  C-,  =  0  (27) 

or 

n  Qil  - 

9i  =  E  7T~  -  C,  =  0  (28) 

1  i=l  Ai  1 

or 

g1  =  { u) t{ S1 }  -  C-,  =  0  (29) 
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dh 
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where  is  the  limiting  value  of  the  displacement. 

Using  Eqs.  26  and  29,  the  Lagrangian  for  a  single-displacement  constrained 
problem  can  be  written  as 


n  /  n  Qil  - 

efAill+x’  ( V ' c’ 


Differentiating  this  equation  with  respect  to  the  design  variable  A^  and 
setting  it  equal  to  zero  gives 

Qil 

Dia1  "  *1  "2  a  0  <31> 

This  equation  can  be  rewritten  as 

2  1 

■  M  ^7  <32> 


1  =  X 


1  „ 

AiPi£i 


A  =  >  _ L! _ 

f  ’  Vt*t 


In  Eq.  35,  e^  is  the  virtual  strain  energy  density.  Eqs.  31  through  35 
represent  the  optimality  criterion  for  a  single-constrained  problem. 

Eq.  35  states  that  in  the  optimum  structure  the  ratio  of  virtual  strain 
energy  density  to  specific  weight  is  the  same  for  all  the  elements. 

The  flexibility  coefficients  are  constant  for  a  statically 
determinate  structure.  In  the  case  of  a  statically  indeterminate  structure. 
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these  coefficients  are  functions  of  the  design  variables.  In  the 
derivation  of  the  optimality  criterion  above,  we  seem  to  have  treated 
Q*1  as  a  constant  with  respect  to  the  design  variable.  However,  even 
without  this  assumption  an  identical  criterion  can  be  derived.  Consider 
the  definition  of  the  displacement  constraint  as  given  in  Eq.  29.  Using 
this  equation  and  Eq.  1,  we  can  write 

C1  =  {ult[K]{s1l  (36) 

where  { s ^ }  is  the  virtual-displacement  vector  associated  with  the  applied 
virtual-load  vector  (S^l.  Differentiating  Eq.  36  with  respect  to  A^  gives 


Eq.  1  can  be  written  as 

<u) t[K]  =  (P)  (38) 


Differentiating  this  equation  with  respect  to  A^  gives 


This  equation  can  be  written  as 


(39) 
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Using  this  relation  and  a  similar  one  for  the  gradient  of  the  virtual  - 
displacement  vector  { s ^  > ^ ,  Eq.  37  can  be  written  as 


3C1 

aX! 


=  -  { 


u!tLj7[K]  {s')  ♦  {u}1 

1  kcW  «*'» 


3A. 


[K] 


is1} 


{ u) 


■  -  {ul 


3  A1 


:k] 


{S1} 


[kL 


Since  [K]  =  '  ,  Eq.  41  becomes 


3C1 

iAT 


1)1 

A1 


(41) 


(42) 


This  derivation  of  the  gradient  of  the  constraint  shows  that  in  deriving  the 
optimality  criterion  it  is  not  necessary  to  assume  as  a  constant. 

Recurrence  Relations 

The  optimality  criterion  derived  in  the  foregoing  section  is  valid 
only  at  the  optimum  and  has  to  be  converted  into  a  recurrence  relation  so 
that  it  can  be  used  in  an  optimization  algorithm.  A  recurrence  relation 
can  be  written  by  multiplying  both  sides  of  Eq.  33  by  A^  and  taking  the 
rth  root.  This  gives 


vi  1 


Aipi*i> 


1/r 

k 


(43) 
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where  k+1  and  k  are  introduced  to  indicate  the  iteration  numbers.  The 
parameter  r  in  Eq.  43  determines  the  step  size.  We  will  call  this  an 
exponential  recurrence  relation.  Eq.  43  can  also  be  written  as 


(44) 


In  this  equation  \-| 


Qi 


i  1 


4>i*i 


is  equal  to  unity  at  the  optimum;  thus. 


/  Qil 

near  the  optimum! A,  - 

v  nVi 


-  l 


is  small  compared  to  unity.  Therefore 


using  the  binomial  theorem  to  expand  the  right  side  of  Eq.  44  and 
retaining  only  the  linear  terms  gives 


(45) 


This  equation  we  will  call  a  linear  recurrence  relation.  In  this  equation 
is  the  error  in  satisfying  the  optimality  criterion. 

The  exponential  and  linear  recurrence  relations  can  be  expressed  also 
in  terms  of  the  virtual  strain  energy  density  e^  used  to  write  the  optimality 
criterion  in  Eq.  35.  To  use  the  recurrence  relation,  one  must  select  the 
proper  step-size  parameter.  For  most  problems,  r=2  is  generally  adequate. 
Increasing  the  value  of  r  reduces  the  step  size. 

The  recurrence  relations  contain  two  unknowns:  the  flexibility 
coefficients  and  the  Lagrange  multiplier  x^ .  The  coefficients 
can  be  determined  by  using  Eq.  10.  The  Lagrange  multiplier  x^  can  be 
obtained  by  using  the  relations  derived  in  the  next  section. 


Ml 


1  2 
Wi 


-  1 
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I 


g 


s 

‘l 


Expressions  for  the  Lagrange  Multiplier 

The  design  variable  Ai  using  Eq.  32  can  be  written  as 


(46) 


The  Lagrange  multiplier 
in  Eq.  28  and  solving  for  x^. 


X.|  can  be  evaluated  by  substituting  Eq.  46 
This  gives 


T  V^ilpi2-i 

i=l 


(47) 


We  can  derive  two  more  expressions  for  the  Lagrange  multiplier  by 
using  the  optimality  criterion  and  the  constraint  equations.  Substituting 
Eq.  34  in  Eq.  26  and  using  Eq.  28,  the  weight  of  the  structure  is  given  by 

W  =  X1C1  (48) 

or 


x 


1 


(49) 


Eq.  48  states  that  at  the  optimum  the  weight  of  the  structure  is  equal 
to  the  product  of  the  Lagrange  multiplier  X^  and  the  limiting  value  of 
the  constrained  . 


Another  expression  for  the  Lagrange  multiplier  can  be  obtained  by 
writing  the  optimality  criterion  as 


«1I 


1  3 

Pi*iA? 


(50) 


4 
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The  Lagrange  multiplier  in  Eqs.  43  and  45  can  be  eliminated  by 
using  one  of  the  expressions  for  A-j  derived  above.  However,  in  the  case 
of  the  single-constrained  problem  if  the  structure  can  be  scaled  by 
using  Eqs.  17  through  22,  a  relative  value  of  the  Lagrange  multiplier  can 
be  used.  Therefore,  one  can  set  A1  equal  to  unity.  The  expressions  for 
the  Lagrange  multiplier  derived  in  this  section  are  strictly  valid  only 
at  the  optimum. 

Effect  of  the  Passive  Elements  on  the  Lagrange  Multiplier 

In  deriving  the  expressions  for  the  Lagrange  multipliers  for  the 
single-displacement  constraint  we  have  assumed  that  all  the  design 
variables  can  be  modified  by  using  the  recurrence  relation.  In  a  prac¬ 
tical  design  problem  there  are  always  constraints  on  the  minimum  and 
maximum  sizes  of  the  elements.  These  elements,  or  other  elements  whose 
sizes  are  not  modified  during  the  iterations  by  using  the  recurrence 
relation  are  called  passive  elements.  One  must  modify  the  algorithm 
to  account  for  these  elements.  When  the  exponential  recurrence  relation 
is  used,  if  the  virtual  strain  energy  associated  with  any  element  is 
negative,  then  that  element  is  sized  by  some  other  criterion  and  becomes 
a  passive  element.  Elements  whose  sizes  are  governed  by  the  recurrence 
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relation  can  be  called  active  elements.  At  the  optimum  the  optimality 
criterion  is  satisfied  only  by  the  active  elements. 

The  contribution  to  the  weight  of  the  structure  and  the  constraint 
can  be  divided  into  two  parts,  one  due  to  the  active  and  the  other  due 
to  the  passive  elements.  Thus  the  weight  of  the  structure  and  the 
constraint  equation  can  be  written  as 


W  =  Z  p  .  A .  £ .  +  W* 
i=l  111 


n1  ^il  _ 

g,  =  Z  t—  +  C*  -  C,  =  0 

1  i=i  Ai  1  1 


where  W*  is  the  weight  of  the  passive  elements  and  is  the  contribution 
of  the  passive  elements  to  the  constraint.  In  Eqs.  53  and  54,  n1  is  the 
number  of  active  elements.  This  change  in  the  definition  of  the  weight 
of  the  structure  and  the  constraint  affects  the  expression  for  the 
Lagrange  multiplier.  With  this  change  Eqs.  47,  49  and  51  can  be  written, 
respectively,  as 


V^1 


f1Vqn0iii 

C,  -  c* 


,  w  -  u 
A1  '  —  * 
C1  -  c, 


c,  -  c; 


1=1  P^A,3 


Hence,  when  there  are  passive  elements  we  should  use  Eqs.  55,  56,  and  57 
instead  of  47,  49,  and  51,  respectively. 
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In  the  case  of  a  single-displacement  constrained  problem  the  explicit 
expression  for  the  Lagrange  multiplier  does  not  have  a  significant 
effect  on  the  algorithm.  However,  in  the  multiple-displacement  constrained 
problem,  the  method  used  to  determine  the  Lagrange  multipliers  plays  a 
major  role  in  the  optimization  algorithm. 

MULTIPLE-DISPLACEMENT  CONSTRAINTS 

The  method  discussed  in  the  last  section  to  design  a  minimum  weight 
structure  would  give  a  minimum  weight  design  only  if  (1)  a  structure  is 
subjected  to  a  single  displacement  constraint,  or  (2)  only  one  constraint 
is  active  at  the  optimum.  The  active  constraints  may  be  defined  as  those 
that  are  satisfied  as  equality  constraints  at  the  optimum.  The  multiple- 
constrained  problem  is  much  more  difficult  to  solve  than  the  single- 
constrained  problem.  The  main  reason  for  this  is  that  a  simple  explicit 
expression  cannot  be  derived  for  the  Lagrange  multipliers.  Also  one  must 
develop  the  logic  to  predict  a  probable  set  of  active  constraints  during 
the  iterations.  If  all  the  active  constraints  are  not  taken  into  con¬ 
sideration,  the  convergence  behavior  is  affected.  If  all  the  constraints 
imposed  on  the  structure  are  considered  potentially  active,  then  the 
algorithm  becomes  inefficient  from  the  point  of  view  of  computational 
effort.  This  is  particularly  true  when  there  are  a  large  number  of 
inequality  constraints  imposed  on  the  structure. 

If  the  structure  is  subjected  to  more  than  one  loading  condition, 
the  loading  conditions  can  be  treated  as  multiple  constraints.  So  to 
simplify  the  notation  we  will  not  consider  the  multiple-loading  conditions 
in  defining  the  problem. 

The  Lagrangian  for  the  multiple-displacement  constrained  problem 
can  be  written  as 


where 


n  m 

L ( A ,  X)  =  z  PiA.*.  +  z  A.g. 

i=l  111  j=l  0  J 


=  Cj 


C.  <  0 

J 


j=l ,  m 


(58) 


(59) 
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n  Qij 
C  =  T. 

J  1=1  N 


C.  =  {ul^S3} 


In  Eqs.  58  and  59,  g.  are  the  inequality  constraints,  x.  are  the 

J  J 

Lagrange  multipliers,  and  m  is  the  number  of  inequality  constraints. 

The  optimality  criterion  can  be  derived  by  differentiating  the  Lagrangian 
with  respect  to  the  design  variables  A^  and  setting  the  resulting  equations 
to  zero.  This  gives 

m  Q.  . 

p  .t  .  -  Z  A  .  I3,  =0  (62) 

j=l  J  A‘ 

where  x.g.=0.  For  the  active  constraints,  x.>0  and  g.=0  and,  for  the 

J  J  J  J 

passive  constraints,  A.=0  and  g./O.  The  optimality  criterion  in  Eq.  62 

\)  J 

can  be  written  in  different  forms  as  follows: 


2  m  Q.  . 

AT  =  z  A.  -4“ 

1  j=l  3  pi*i 


m  Q.  . 

A.  =  Z  X . 

1  j=l  J  pi£i Ai 


m  Q. 


1  =  Z  X 


i  2 
j*i  °ic'iAi 


m  e.  . 
1  ■  z  a. 
j=l  3  pi 
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or 


or 


where 


.  Ml 

{u}J[k]i  z  k  .{sJ}. 


1  = 


ill 


piS‘Ai 


£ 

J=I 


A,{S3}. 


Ai£i 


(67) 


(68) 


(69) 


The  optimality  criterion  as  written  in  Eqs.  66  and  68  has  a  physical 

interpretation.  According  to  Eq.  66,  at  the  optimum  the  weighted  sum  of 

the  ratio  of  virtual  strain  energy  density  to  mass  density  is  equal  to 

unity  for  all  the  elements,  where  the  weighting  parameters  are  the 

Lagrange  multipliers.  The  virtual  strain  energy  density  in  each  element 

is  due  to  a  virtual  load  vector  associated  with  each  distinct  constraint. 

Eq.  68  which  is  also  an  optimality  criterion  states  that  the  virtual 

strain  energy  to  mass  density  is  equal  to  unity  for  all  elements  where 

m 

a  single  virtual  load  vector  equal  to  £  x  .{ SJ }  is  applied  to  the 

j=l  J 

structure.  In  using  the  recurrence  relation  based  on  Eq.  65  or  66,  one 
must  determine  in  each  element  the  virtual  strain  energy  due  to  m  virtual 
loads.  But  if  the  recurrence  relation  based  on  Eq.  67  is  used,  it  is 
only  necessary  to  consider  one  virtual  load  vector  to  determine  the  virtual 
strain  energy  in  an  element.  However,  one  must  know  the  actual  or  relative 
values  of  the  Lagrange  multipliers  before  a  single  virtual  load  vector 
can  be  assembled.  This  is  not  always  possible  as  we  will  see  when  we 
discuss  the  methods  to  determine  the  Lagrange  multipliers. 

The  ' n *  optimality  conditions  and  the  'm'  constraints  must  be 
satisfied  by  the  optimum  design.  These  are  (m+n)  nonlinear  equations 
that  must  be  solved  to  determine  the  'n1  design  variables  and  the  'm' 
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Lagrange  multipliers.  The  optimization  algorithm  solves  these  nonlinear 
equations  by  an  iterative  scheme.  In  an  iterative  method  the  design 
variables  are  modified  by  using  the  recurrence  relation  so  that  the 
optimality  criterion  is  satisfied.  The  Lagrange  multipliers  are  evaluated 
on  the  basis  that  the  constraint  conditions  are  satisfied  when  the  design 
variables  are  changed.  Since  the  Lagrange  multipliers  and  the  design 
variables  are  functions  of  each  other,  any  change  in  one  nonlinearly 
affects  the  other,  and  this  necessitates  the  use  of  an  iterative  method. 


Recurrence  Relations  for  Multiple  Constraints 

The  recurrence  relations  for  the  multiple  displacement  constrained 
problem  can  also  be  divided  into  two  categories:  (1)  an  exponential 
form,  and  (2)  a  linear  form,  as  we  have  done  for  the  single-displacement 
constrained  problem.  The  method  of  deriving  them  is  similar  to  the  one 
we  used  to  write  Eqs.  43  and  45  for  the  single-constrained  problem. 

An  exponential  recurrence  relation  can  be  written  by  multiplying 
Eq.  65  by  A^  and  taking  the  rth  root.  This  gives 


A  . 
0 


*LL 


1/r 


pi*iAi 


(70) 


where  k+1  and  k  are  the  iteration  numbers.  In  Eq.  70  the  parameter  r 
▼  > 

determines  the  step  size.  Using  the  same  argument  given  to  derive 
Eq.  45,  we  can  write  a  linear  recurrence  relation  for  the  multiple- 
displacement-constrained  problem.  This  gives 


A, 


j  »i'X 


This  equation  can  also  be  written  as 


Ak+1  =  Ak  +  1 

«i  +  r 


A, 


j  P/iAi 


k 


(71) 


(72) 
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In  this  equation  the  term  in  parenthesis  is  the  correction  needed  to  the 
design  variable  to  satisfy  the  optimality  criterion  as  defined  in 
Eq.  64.  Thus  we  can  say  that  the  objective  of  modifying  the  design 
variable  by  using  the  recurrence  relation  is  to  move  the  current  design 
towards  a  design  that  satisfies  the  optimality  criterion.  When  the 
step-size  parameter  r=l ,  Eqs.  70  and  71  reduce  to  the  same  equation. 

For  the  multiple-constrained  problem,  as  in  the  case  of  a  single  con¬ 
straint,  r= 2  can  be  considered  as  the  normal  step  size,  but  some  problems 
require  one  to  increase  r  to  reduce  the  step  size.  In  the  exponential 
recurrence  relation  the  design  variable  is  modified  by  multiplying  it 
by  a  quantity  which  is  equal  to  unity  at  the  optimum,  and  in  the  linear 
recurrence  relation  the  design  variable  is  modified  by  adding  a  quantity 
which  is  equal  to  zero  at  the  optimum.  Eqs.  70  and  71  can  also  be  written 
by  using  the  optimality  criterion  defined  by  Eqs.  66  and  68. 

Methods  to  Determine  the  Lagrange  Multipliers 

In  order  to  use  any  one  of  the  recurrence  relations  the  coefficients 
and  the  Lagrange  multipliers  have  to  be  known.  The  coefficients 
Qi j  can  be  determined  by  using  Eq.  10.  The  equations  to  determine  the 
Lagrange  multipliers  can  be  derived  by  considering  the  effect  of  a  change 
in  the  design  variables  or  Lagrange  multipliers  on  the  constraint  equations. 
We  will  derive  in  this  section  different  methods,  some  rigorous  and  some 
not  rigorous.  The  rigorous  methods  give  equations  that  are  reliable  but 
the  methods  need  more  computational  effort.  The  nonrigorous  methods  are 
simple  and  need  less  computational  effort,  butgive  algorithms  which  may 
not  lead  to  a  minimum  or  are  slow  to  converge.  In  discussing  the  different 
methods  we  will  point  out  their  advantages  and  disadvantages  and  their 
relationship  to  one  another. 

Recurrence  Relations  for  the  Lagrange  Multipliers 

A  recurrence  relation  to  estimate  the  Lagrange  multipliers  can 
be  written  by  assuming  that  all  the  constraints  in  Eq.  59  are  equality 
constraints.  This  gives 

Cj  -  tj  (73) 
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Multiplying  both  sides  of  this  relation  by  x^  and  taking  the  bth  root,  a 

3 

recurrence  relation  can  be  written  as 


where  k  refers  to  the  iteration  number  and  the  parameter  b  controls  the 
step  size.  The  advantages  of  using  this  recurrence  relation  are: 

(1)  An  equivalent  single  virtual  load  vector  can  be  used  to 
determine  the  virtual  strain  energy  in  an  element,  and  individual  values 
of  Q.  .  for  different  constraints  need  not  be  determined. 

*  J 

(2)  It  is  not  necessary  to  predict  which  constraints  are 
potentially  active.  Repeated  use  of  Eq.  74  reduces  the  value  of  a 
Lagrange  multiplier  corresponding  to  a  passive  constraint,  and  its 
contribution  to  the  virtual  load  vector  is  reduced  after  each  iteration. 


(3)  The  computational  effort  required  to  determine  the 
Lagrange  multipliers  and  the  virtual  strain  energy  in  an  element  is 
minimal  as  compared  to  other  methods. 

The  disadvantages  are  (1)  Initial  values  of  the  Lagrange 
multipliers  have  to  be  assumed,  and  (2)  Convergence  to  the  minimum-weight 
design  is  slow  and  large  oscillations  sometimes  occur  in  the  scaled 
weight  of  the  structure. 


The  recurrence  relation  in  Eq.  74  can  be  written  as 


C. 


Since  at  the  optimum  ^  is  nearly  equal  to  unity,  the  difference^ 

is  small  as  compared  to  unity.  We  can  therefore  expand  Eq.  75  by' using' 
the  Binomial  theorem  and  retain  only  the  linear  terms.  This  gives 


■(I: ) 


k+1  -  xk 

J  J 


(76) 
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This  equation  can  also  be  written  as 


We  will  show  in  the  next  section  that  this  linear  recurrence  relation 
for  the  Lagrange  multipliers  is  an  approximation  to  a  set  of  linear 
equations  that  can  be  used  to  determine  the  Lagrange  multipliers. 

Linear  Equations  to  Determine  the  Lagrange  Multipliers 
A  set  of  equations  to  determine  the  Lagrange  multipliers  can  be 
obtained  by  considering  the  change  in  the  constraint  due  to  a  change  in 
the  design  variable  A..  The  change  in  the  jth  constraint  can  be  written 

J 

as 

Ag^  =  gj (A  +  aA)  -  g^(A)  (78) 


n  3g . 

■  £  (79) 

1=1  3fli  ' 


If  the  change  in  the  design  variable  is  assumed  to  be  such  that  the 
constraints  at  point  { A+dA}  are  satisfied  i.e.  g.(A+dA)=0,  then  using 

J 

Eq.  59,  Eq.  79  can  be  written  as 


n 

£ 

i  =1 


(80) 


where  k  refers  to  the  iteration  number.  Differentiating  Eq.  60  with 
respect  to  the  design  variable  A-  and  assuming  that  is  a  constant 
gives 


3g . 
_1L 
3A . 


(81) 


The  change  in  the  design  variable  from  the  kth  to  the  (k+l)th  iteration 
is  given  by 

AAk  =  Ak+1  -  Ak  (82) 

ill 
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Substituting  A./  from  Eq.  71  in  this  equation  gives 


/  “  J  Q<i  .  A  Ak 


Using  Eqs.  81  and  83,  Eq.  80  can  be  written  as 


;  H ^ 

3=1  p  \  i-i  p,*X 


k  -  n  /Qi  i 

■ r  (cj  •  *  dt 


Since 


,"=,  (^)k  *  CJ' E<" 


84  can  be  written  as 


£  \k+]  (  £  QljQl'§  V  (r+l)Ck  -  rC 
P=i  p  Vi=i  pi£iA?y  J 


These  are  'm'  equations  corresponding  to  the  'm'  displacement  constraints, 
which  can  be  used  to  determine  the  Lagrange  multipliers  in  the  iterative 
algorithm.  At  each  iteration  one  must  use  only  those  equations  corre¬ 
sponding  to  the  active  constraints,  giving  positive  Lagrange  multipliers. 
This  can  generally  be  achieved  by  considering  only  those  constraints 
that  are  closest  to  the  constraint  surface  and  by  eliminating  the 
equations  yielding  negative  Lagrange  multipliers. 

When  there  are  passive  elements  whose  sizes  are  not  governed  by  the 
optimality  criterion,  we  have  to  separate  the  contribution  of  those  ele¬ 
ments  in  the  summation  of  Eq.  84.  This  gives 


i  \ 
p=l 


k+l/  Jl  Qi,jQip\  =  {ck  .  c  )  +  "if 

r<Cj  Cjl  44 


n  /q..  .  \  k 

-  r  £  f  )  A\(p) 

’"V'\Ai  / 
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k  k  k 

where  'n^'  is  the  number  of  active  elements  and  AA^p)  =  A^p)  -  Ai . 

1/ 

A ^ ( p )  is  the  size  of  a  passive  element  dictated  by  considerations  other 
than  the  optimality  criterion.  Generally  the  passive  elements  are  those 
whose  sizes  are  governed  by  a  minimum  or  a  maximum  size.  For  a  problem 
with  stress  and  displacement  constraints,  if  we  decide  to  treat  the  stress 
constraints  as  passive  constraints,  then  the  elements  whose  sizes  are 
governed  by  stresses  will  be  included  in  the  passive  category. 

The  advantages  of  using  Eqs.  85  and  86  to  determine  the  Lagrange 
multipliers  are: 

(1)  It  is  not  required  to  assume  initial  values  of  the  Lagrange 
multipliers,  since  the  Lagrange  multipliers  are  evaluated  by  solving  the 
equations . 

(2)  Since  the  equations  contain  coupling  terms  which  take  into 
account  the  interdependence  of  the  different  constraints,  the  Lagrange 
multipliers  are  more  realistic.  This  is  particularly  true  in  the  case 
of  a  structure  where  the  constraints  are  sensitive  to  design  changes. 

(3)  The  convergence  behavior  is  generally  better  than  other 
approaches . 

The  disadvantages  are: 

(1)  A  single  virtual  load  vector  cannot  be  assembled  since  the 
Lagrange  multipliers  are  not  known  before  the  flexibility  coefficients 
Q.  .  corresponding  to  the  constraints  are  determined. 

'  J 

(2)  The  computational  effort,  to  determine  the  flexibility 
coefficients  and  to  assemble  the  coefficients  of  the  Lagrange  multipliers 
in  Eqs.  84,  85,  and  86  substantially  increases  with  an  increase  in  the 
number  of  potentially  active  constraints. 

(3)  It  is  necessary  to  solve  a  set  of  simultaneous  equations. 

(4)  Some  scheme  has  to  be  used  in  order  to  determine  the  potentially 
active  constraints  and  to  eliminate  those  equations  corresponding  to  the 
passive  constraints. 
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At  the  optimum  the  active  constraints  are  satisfied  as  equality 
constraints,  i.e.  C.=C.  and  Eq.  85  becomes 

J  J 


p  V-'  “A*? 


This  relation,  although  valid  only  at  the  optimum,  can  be  used  to 

estimate  the  Lagrange  multipliers.  Eq.  87  can  also  be  derived  directly 

by  using  the  optimality  criterion  and  the  constraint  equations.  The 

c 

optimality  criterion  can  be  written  as 

{1}n  =  ^nxm  U}m  (88) 

The  constraint  equations  can  be  written  as 

[F]mvn  {lK  =  (C)  (89) 

mxn  n  m 

In  Eqs .  88  and  89,  {1}  is  a  vector  with  all  the  elements  equal  to  unity. 
Using  these  equations  we  can  write 

[F]  [Q]  UK  =  (Cl  (90) 

mxn  ^Jnxm  m  m 


tRU  “>■  ’  <91 

where  [R]mxm  -  ^^mxn  ^nxm 

In  Eq.  91  the  elements  of  [R]  are  given  by 

n  Q. .Q. 

Rin  “  2  (* 

JP  i-i 

Comparing  Eqs.  87  and  91  it  is  seen  that  they  are  identical. 

Eq.  87  can  also  be  derived  by  writing  the  optimality  criterion  as 


L  A.  3 
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and  substituting  it  into  the  constraint  equation 

n  Qij  - 
I  =  C  . 


A 

i  =1  Ai 


The  recurrence  relation  for  the  Lagrange  multipliers  given  by 
Eq.  76  can  be  shown  to  be  an  approximation  to  the  linear  equations  in 
Eq.  85.  The  matrix  multiplying  the  vector  {X}  in  Eq.  85  is  square. 

If  we  neglect  the  off-diagonal  terms,  these  equations  become 

Ak+1  [!  =  ( r+1  )C^  -  rC .  (< 

J  i  =1 V  P..2...A;/  J  J 


These  equations  are  uncoupled  and  assume  that  the  constraints  are  inde¬ 
pendent  of  one  another.  With  this  assumption,  and  using  the  optimality 

o  a 

criterion,  (Q.  ./p. 8. - A? ).  in  Eq.  95  can  be  replaced  by  1/x^.  Also 

I J  I  I  1  K  J 


recal ling  that 


n  Qij\  k 

l  a-*1)  =  c£,  Eq.  95  can  be  written  as 

Ai  Jk 

xk+1 

-V  ci =  (r+i)ci  -  rCi 

A.  J  J  J 

J 


This  equation  would  be  identical  with  Eq.  76  if 

l-r  Of 

This  shows  that  Eq.  76,  which  is  a  linearized  form  of  Eq.  74  is  an 
approximation  to  Eq.  85. 

Newton-Raphson  Method 

An  iterative  algorithm  to  solve  the  optimality  criterion  and 
the  constraint  equations  can  be  developed  by  considering  the  change  in 
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a  constraint  due  to  a  change  in  the  Lagrange  multipliers.  The  change 
in  a  constraint  can  be  written  as 


m  3g . 


Ag.  =  g.  (a  +  &A)  -  g . (x )  =  e  -rr* 
J  J  J  j=l  dA 


AA 


(99) 


Since  in  the  Newton  Raphson  method  the  change  aa  is  selected  to  satisfy 
the  condition  g.(A+AA)=0,  using  Eq.  99  an  iterative  relation  can  be 

J 

written  as 


-  <t>{ g> k  =  [H]{ Ak+1  -  Ak} 


(100) 


or 


{A}k+1  =  ( A } k  -  ♦  [H]“1{g}|c 


(101) 


where  [H]”^  is  the  inverse  of  the  Hessian  [H]  whose  elements  are 

39,- 


Hjp  =  3X 


(102) 


and  i?)  is  a  parameter  introduced  to  control  the  step  size. 


Differentiating  Eq.  60  with  respect  to  A.,  recalling  that  A.  is 

J  J 

related  to  A.  by  Eq.  63  and  assuming  Q.  .  as  constants,  we  can  write 

!  I  J 


i ; 

2  i-l 


Q.  .0. 

’1*1*1 


\  * 


93) 


Eq.  101  can  be  used  to  update  the  initially  assumed  Lagrange  multipliers. 
The  summation  in  Eq.  103  is  carried  out  only  over  the  active  elements. 

The  iterative  process  in  this  method  consists  of  using  Eq.  101  to  esti¬ 
mate  the  Lagrange  multipliers  and  Eq.  63  to  modify  the  design  variables 
alternatively  until  the  constraint  equations  are  satisfied  within  a 
specified  tolerance.  The  disadvantages  of  this  method  are:  (1)  It  is 
essential  to  assume  initial  values  of  the  Lagrange  multipliers,  (2)  The 
passive  constraints  cannot  be  easily  separated  from  the  active  constraints. 
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and  (3)  If  too  many  iterations  are  performed  to  satisfy  the  constraints, 
the  design  might  move  away  from  the  region  where  the  coefficient 
are  valid. 


The  Newton-Raphson  algorithm,  can  be  shown  to  be  related  to  the 
linear  equations  in  Eq.  85.  Eq.  101  can  be  written  as 


[HJk{X}k+1  =  [H]k{A}k  -  <|>{ g>k 


(104) 


Using  Eq.  103,  this  equation  can  be  written  as 


m  .  .  n  / Q.  .0. 

E  Xk+'  e/-^ 

P-1  P 


m 

e  x; 


p-l  p1=llp  ,*,*?)  J 


(105) 


Substituting  Af  from  Eq.  63  into  Eq.  60  and  rearranging,  one  obtains 


gk  =  J  Ak  J  /Q1JQ1P 

9j  P=iXpiMPiv? 


-  c 


(106) 


b  _ 

Using  this  equation  and  recalling  that  g.  =  (C.  -  C.),  Eq.  105  can  be 

J  J  J 

written  as 


m  k+1  n 

E  x”  E 

P=1  P 


i=l 


=  (2<j>+l  )Ck  -  2*Cj 


(107) 


This  equation  is  identical  to  Eq.  85,  for  2<p=r. 


The  advantage  of  using  Eq.  107  in  the  Newton-Raphson  method  instead 
of  Eq.  101  is  that  one  need  not  assume  the  initial  values  of  the  Lagrange 
multipliers.  Eq.  107  can  be  solved  at  each  iteration.  This  approach  we 
re*er  to  as  the  modified  Newton-Raphson  method. 
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SECTION  IV 
STRESS  CONSTRAINTS 


An  algorithm  can  be  developed  to  design  a  structure  with  stress 
constraints  by  using  two  approaches.  The  first  approach  is  to  specify 
the  constraint  as  a  function  of  the  stress  in  an  element,  and  the  second 
is  to  convert  the  stress  constraint  into  an  equivalent  displacement  con¬ 
straint.  These  two  approaches,  depending  on  the  approximations  used  to 
derive  the  criterion,  lead  to  different  algorithms.  For  convenience  we 
will  consider  only  bar  elements,  but  the  conclusions  and  the  algorithms 
can  be  easily  extended  to  other  membrane- type  elements. 


FULLY  STRESSED  DESIGN  (FSD)  METHOD 

The  stress  constraint  in  the  ith  element  can  be  written  as 


(108) 


where  is  the  actual  stress  and  is  the  maximum  allowable  stress  in 
the  ith  element.  The  number  of  constraints  is  equal  to  the  number  of 
elements.  The  stress  in  the  ith  element  is  given  by 


(109) 


where  T^  is  the  force  in  the  ith  element.  Using  Eqs.  26,  108,  and  109, 
the  Lagrangian  for  the  stress-constrained  problem  can  be  written  as 


n  n  /  T. 

L(A,  x)  =  z  p • A. i •  +  z  X.[  — 

1=1  1=1  \Vi 

Differentiating  this  equation  with  respect  to  the  design  variable  A^ 
and  setting  it  equal  to  zero  gives 


pi*i 


x. 


i 


Aft 


Vi 


n 

z 

p=l 


x  ^ 
XP  3A. 


(Ill) 
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3TP 

In  this  equation  the  term  ,  which  is  the  gradient  of  the  force  in  a 

bar,  cannot  be  explicitly  written.  For  a  determinate  structure,  this 
gradient  is  zero,  since  the  force  in  a  bar  is  independent  of  the  areas  of 

n  3T 

the  areas  of  the  elements.  If  we  assume  that  z  X  =  0,  which  is 

P=i  "  “Ai 

not  true  for  an  indeterminate  structure,  Eq.  Ill  becomes 

T. 

pi*i  ‘  xi  =  0  012) 

A,  a. 


or 


Now  if  the  stress  constraints  imposed  on  the  structure  are  assumed  to  be 
satisfied  as  equality  constraints,  then  using  Eqs.  108  and  109  we  can 
write 


(114) 


Substituting  Eq.  113,  in  this  equation  and  solving  for  gives 


(115) 


Substituting  this  relation  in  Eq.  113  gives 


(116) 


This  equation  states  that  at  the  optimum  the  area  of  the  element  is 
equal  to  the  force  in  the  bar  divided  by  the  maximum  allowable  stress  in 
the  bar,  or  the  stress  in  an  element  is  equal  to  the  maximum  allowable 
stress  in  that  element.  This  is  the  well-known  Fully  Stressed  Design 
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optimality  criterion, 
can  be  written  as 


A  recurrence  relation  based  on  this  criterion 


At  the  optimum  T\  =  A^.  Substituting  this  in  Eq.  113  we  can  write 


x,  -  A1Pfei 

=  weight  of  the  element 


(118) 


The  use  of  the  Fully  Stressed  Design  algorithm,  if  applied  to  an  inde¬ 
terminate  structure,  is  an  approximation,  since  for  this  structure 

n  3T 

l  A  — ^  o.  Therefore  the  FSD  algorithm  does  not  necessarily  give 
p-1  p  3Ai 

a  minimum  weight  design  for  an  indeterminate  structure,  and  this  is 
particularly  true  when  the  allowable  stress  in  all  the  elements  is  not 
the  same.  The  design  obtained  by  the  Fully  Stressed  Design  algorithm 
for  an  indeterminate  structure  with  unequal  stresses  sometimes  not  only 
gives  a  non-optimum  design  but  also  gives  a  design  with  a  bad  distribution 
of  material  and  an  inefficient  load  path. 


STRESS  CONSTRAINT  WITH  EQUIVALENT  DISPLACEMENT  CONSTRAINT 
The  stress  in  the  ith  bar  element  can  be  written  as 

°i  "  [0n  d123i  ru>l  Ci?) 

u2 

L.  _lj 


1  2 

where  U  and  U  are  the  longitudinal  displacement  of  the  two  nodes 

E.  E. 

defining  the  ith  bar.  In  Eq.  119,  D^  =  -  •—  and  D^  =  where  Ei  is 

the  Young's  modulus  and  is  the  length  of  the  bar.  The  stress  con¬ 
straint  in  the  ith  element  can  now  be  written  as 


g.  =  o1  -  o1  <  0 

(120) 

or 

=  (D11U1  +  D12U2).-  ^<0 

(121) 
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Thus  we  have  replaced  the  stress  constraint  in  an  element  by  an  equivalent 
constraint  on  the  linear  combination  of  the  displacements  in  the  longi¬ 
tudinal  direction  at  the  two  nodes  defining  the  element.  The  constraint 
equation  can  now  be  written 


or 


or 


where 


g .  =  a  .  -  a  .  <  0 
J  J  ~ 


j=l  . .n 


n  R.  . 

q  .  =  i  -  a  .  <  0 

J  i=i  Ai  J  - 


g,  =  {u}t{SJ)  -0,10 

J  J 


Rij  ' 


(122) 

(123) 

(124) 

(125) 


There  will  be  n  constraints  equal  to  the  number  of  elements  in  a  bar 
structure.  In  Eqs.  125  {sJ)1-  is  the  vi rtual -di spl acement  vector  associated 
with  the  ith  element  due  to  the  virtual  load  vector  (SJ),  which  is 
equivalent  to  the  forces  D,,  (=  -  E./£.)  and  D,2  (=E./£.)  acting  in  the 
longitudinal  direction  at  the  two  nodes  at  the  ends  of  jth  element. 

Using  Eqs.  26  and  123  the  Lagrangian  for  the  stress-constrained  problem 
can  be  written  as 


L(A,  X) 


n 

I  p.A.£. 

1-1  1  1  1 


n 


+  z 
j=l 


(126) 


Differentiating  this  relation  with  respect  to  the  design  variable  and 
setting  the  result  equal  to  zero  gives 


p/i 


n 

z 

j=l 


X  .  g 


R.  . 
_LL  = 


or 


(127) 


1 


*  2 
j=l  J  p,£,.Af 


R.  . 

1L 


i  ^ 


(128) 
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where 


x  .  >_  0  and  x  .a  .  =  0 
J  J  J 


(129) 


This  optimality  criterion  derived  for  the  stress  constraints  has 
the  same  form  as  the  optimality  criterion  in  Eq.  65  for  the  multiple- 
displacement-constrained  problem.  Hence  following  the  same  procedure 
and  logic  we  can  derive  the  recurrence  relations  and  the  equations  to 
determine  the  Lagrange  multipliers.  The  recurrence  relations  for  the 
stress-constrained  problem  equivalent  to  Eqs.  70  and  71  are 


»r  -  <  ( "  ^ 

v=1  J 


(130) 


Ai+1  =  Ai  I1  +  r 


/  n  R.  .  \ 

(  s  A 

\j=l  J  P^Af  / 


(131) 


Similarly  the  equations  to  determine  the  Lagrange  multipliers  equivalent 
to  Eqs.  74,  76,  77,  85,  and  87,  respectively,  are 


xk+1  =xk 

J  J \  b  bo. 


/  xk+1 

k  .  k  /.  j 

a  ■  -  o  .  -  ba  .  I  1  -  , 

3  \  Xj 


(132) 


(133) 


(134) 


n  /R.  .R.  \ 
z  (JLJ  IP.  \ 

1=1VlVi7l 


=  (r+1 )o j  -  roj 


(135) 
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In  the  case  of  a  stress-constrained  problem,  the  number  of  constraints 

is  equal  to  the  total  number  of  stress  constraints  imposed  on  the 

structure.  For  a  bar  structure  the  number  of  constraints  will  be  equal 

to  the  number  of  elements.  If  the  structure  is  idealized  with  membrane 

elements  and  stress  constraints  are  imposed  on  a  ,  a  ,  a  ,  there  will  be 

x  y  xy 

three  constraints  associated  with  each  element. 

In  Eq.  137,  if  the  off-diagonal  terms  multiplying  the  vector  of 
Lagrange  multipliers  are  neglected,  this  equation  can  be  approximated  as 


x 


k+1 

P 


n  /R.  R.  \ 

"Vi'iA7k 


(138) 


In  addition  if  we  assume  that  the  virtual  strain  energy  in  an  element  is 
only  due  to  the  virtual  load  in  that  element,  which  is  true  for  a 
statically  determinate  structure,  then  Eq.  138  becomes 


x 


k+1 

i 


(139) 


Since  in  this  equation  we  are  assuming  that  the  elements  are  not 
T.n. 

interdependent,  1  .  In  addition  if  we  assume  that  the  elements 
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T. 

are  fully  stressed,  i.e.  k.  =  ,  a  simple  expression  for  the  Lagrange 

°i 

multiplier  can  be  written  as 


x 


k+1 

i 


(140) 


We  found  that  this  simple  expression  leads  to  a  correct  near-minimum- 
weight  design  for  the  case  of  different  allowable  stresses  in  a  structure, 
if  the  stresses  in  the  elements  at  the  optimum  are  equal  to  the  allowable 
stresses  in  the  elements. 


The  Newton-Raphson  algorithm  derived  for  the  multiple-displacement 
constrained  problem  can  also  be  derived  for  the  stress-constrained 
problem  by  using  the  same  procedure.  The  iterative  relation  would  be 

Ulk+1  =  <X>*  -  ♦CH]-1{g}k  (141) 

where 


i 


i 


S’ 

i=l 


Vie 


(142) 


The  optimality  criterion  for  the  stress  constraints,  Eq.  128,  can 
also  be  written  as 


or 


,  .  "  VulS[k¥sJ)i 

j-1  Vi'i 


{u)f[kL  "  A  (Sj}. 

1  j=l  J  1 

PiVi 


(143) 


(144) 


i 

In  Eq.  144,  if  the  resultant  displacements  E  x.{sJ).  associated  with 

j=l  J 

the  ith  element  due  to  the  virtual-load  system  are  equal  to  {u)i  owing 
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to  the  applied  load,  the  optimality  criterion  becomes 


{u}J[k].{u}. 

Pi£iAi 


(H5) 


n  i 

The  assumption  that  I  \.{s  }.  =  {ul.  is  equivalent  to  assuming 

j=l  J  1 

that  the  product  of  the  Lagrange  multiplier  and  the  virtual  load  applied 
to  each  element  is  equal  to  the  force  in  the  element  owing  to  the  load 
{ P} .  Eq.  145  can  also  be  written  as 


«iVi 


(146) 


where  is  the  strain  energy  stored  in  the  element  and  e..  is  the  strain 
energy  density. 

This  criterion  states  that  at  the  optimum  the  ratio  of  the  strain  energy 
density  to  the  mass  density  is  the  same  for  all  the  elements.  The 
criterion  corresponds  to  a  constraint  on  the  generalized  stiffness  and 
can  be  derived  by  replacing  the  virtual  load  vector  {S1}  by  the  actual 
load  vector  (Pi  in  Eq.  29. 

A  recurrence  relation  based  on  Eq.  146  can  be  written  as 


A*+1  =  Aj 


(147) 


This  recurrence  relation  can  be  used  as  an  approximation  to  Eq.  130  to 
solve  a  stress-constrained  problem.  In  this  case  also  if  the  maximum 
allowable  stress  in  the  different  elements  is  not  the  same,  it  may  not 
lead  to  a  minimum-weight  design.  The  advantage  of  Eq.  147  is  that  it 
is  not  required  to  calculate  the  virtual  strain  energies  in  the  elements 
and  the  Lagrange  multipliers.  A  linear  form  of  Eq.  147  can  be  written  as 


(148) 
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SECTION  V 

ALGORITHM  WITH  THE  RECIPROCAL  DESIGN  VARIABLE 


In  Sections  3  and  4  we  have  used  the  direct  design  variable  A^  to 
define  the  objective  function  and  the  constraint  equations.  The 
recurrence  relations  and  the  equations  to  determine  the  Lagrange 
multipliers  were  expressed  in  terms  of  A^ .  Now  we  will  see  how  the 
optimality  criterion  and  the  algorithm  is  affected  by  defining  the 
problem  in  terms  of  the  reciprocal  design  variable  z^ .  The  reciprocal 
design  variable  z^  and  the  direct  design  variable  A^  are  related  by 
Eq.  9.  We  will  derive  the  equations  only  for  the  multiple-displacement- 
constrained  problem,  since  the  conclusions  from  this  can  be  readily 
extended  to  the  stress-constrained  problem. 


OPTIMALITY  CRITERION 

The  weight  of  the  structure  can  be  written  as 


n  P,*, 

W(z)  =  E  4-1 


i  =1  i 


The  constraint  equations  in  terms  of  zi  can  be  written  as 

j=l .m 


g .  =  C .  -  C .  <  0 
J  J  J  ~ 


where 


C .  =  e  Q.  .z. 

J  i=l  1J  1 

Using  Eqs.  26  and  150,  the  Lagrangian  can  be  written  as 
L(z,  A) 


n  p  £  m  /  m  _\ 

^  4?,  v.  -  h) 


(149) 


(150) 


(151) 


(152) 


Differentiating  this  equation  with  respect  to  the  design  variable  z^  and 
setting  the  result  equal  to  zero  gives  the  optimality  criterion, 
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I 


I 


I 


I 


or 


1  = 


Vi 


m 


L  A.Q..Z. 


j  =  l 


J  1J  1 


where  A  .>0  and  A.g.=0. 

3  3  3 

Eq.  154  can  also  be  written  as 


1 


m 

T.  A  .e.  . 


j=l 


J  iJ 


054) 


(155) 


where  e^  .  is  the  virtual  strain  energy  density.  Comparing  the  optimality 
criterion  for  the  direct  design  variable  in  Eq.  66  with  Eq.  155  shows 
them  to  be  equivalent,  even  though  one  is  the  reciprocal  of  the  other. 


RECURRENCE  RELATIONS 


Applying  the  same  logic  used  to  write  the  exponential  recurrence 
relation  for  Ai ,  (Eq.  70),  the  recurrence  relation  for  zi  can  be  written 
by  using  the  optimality  criterion  in  Eq.  154.  This  gives 


2 

.  ,  /  m  A  .Q.  .z .  \ 

.  2k  (  J  j!lli  ) 

1  '  \j=l  °1l1  / 


This  recurrence  relation  is  equivalent  to  Eq. 
obtained  from  the  other  by  using  Eq.  9. 


-1/r 

k 

70  since  one  can  be 


(156) 


The  linear  recurrence  relation  for  the  reciprocal  design  variable 
can  be  obtained  by  expanding  Eq.  156  using  the  binomial  theorem  and 
retaining  only  the  linear  terms.  This  gives 


2 

A.Q.  .zT 
J  i 

pi*i 


(157) 


This  equation  can  be  expressed  in  terms  of  A^  by  using  Eq.  9.  This  gives 


(158) 
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Comparing  this  equation  with  the  linear  recurrence  relation  for  A.  in 
Eq.  71  shows  that  they  are  not  equivalent. 

EVALUATION  OF  THE  LAGRANGE  MULTIPLIERS 

The  set  of  equations  to  determine  the  Lagrange  multipliers  can  be 
obtained  by  considering  a  change  in  the  constraint  owing  to  a  change  in 
the  design  variable  z^ .  The  change  in  a  constraint  can  be  written  as 

Agj  =  gj  <z  +  Az)  _  9j (159) 
n  ag 

=  T  — J-  az.  (160) 

i=l  azi  1 


Using  Eqs.  150  and  157  and  remembering  that  g.(z+Az)=0,  Eq.  160  reduces  to 


/  m  Qi  iQi  Dzi  \  k  -  n 

~^r)  =  r(cj  •  v  *  4,  (Vi>i 


(161) 


Comparing  this  relation  with  Eq.  84  shows  that  they  are  equivalent  and 
that  one  can  be  obtained  from  the  other  by  using  the  relationship  between 
Aj  and  z^ .  Similarly  one  can  show  that  the  iterative  relations  for  the 
Newton-Raphson  method  for  z.  are  also  equivalent  to  those  of  the 
direct-design  variable  A^  (Eq.  101). 


The  equations  equivalent  to  Eq.  156  and  157  for  the  stress-constrained 
problem  can  be  written  as 


(162) 


The  equations  derived  in  this  section  show  that  the  definition  of  the 
problem  in  terms  of  the  reciprocal  design  variable  affects  the  linear 
recurrence  relation.  However,  the  optimality  criterion,  the  exponential 
recurrence  relation  and  the  equation  to  determine  the  Lagrange  multipliers 
are  not  changed  so  as  to  affect  the  behavior  of  the  algorithm. 
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SECTION  VI 

SYSTEM  STABILITY  CONSTRAINT 

The  constraint  equation  for  the  linear  static  buckling  of  a 
structure  can  be  written  as 

90  ‘  vi  -  » 


(164) 


where  y  is  the  lowest  critical  load  factor  and  y.  is  given  by  Eq.  16. 

J 

The  Lagrangian  for  the  constraint  on  static  stability  can  be  written  as 


L(A,  x)  =  £  P-iM,-  +  E  -  u) 

i=l  111  j=l  J  J 


(165) 


The  optimality  criterion  is  obtained  by  differentiating  Eq.  165  with 
respect  to  Ai  and  setting  the  result  equal  to  zero.  This  gives 


m  3y  . 

p  .  i .  +  Z  \  .  =  0 

The  gradient  of  y.  can  be  found  by  writing  Eq.  15  as 
J 

+  u -{n}^[Kr]{n}  .  =  0 
J  /  ■  J  J  J  6  J 


(166) 


(167) 


and  differentiating  it  with  respect  to  A^ .  This  gives 

^  if  It  t  a 

Slf  ■  •  A7  (nJjtkljCnlj  -  UjtnJj  5J7 


S557[KG]  <n,J 


(168) 


The  second  and  third  term  on  the  right  side  of  this  equation  are  equal 
to  zero,  and  therefore  we  can  write 


^ . .  j_  (nlk[l°ilnl.i 

A.  Ai  {n}^[Kr]{n}  • 
J  6  J 


(169) 
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w. 


N 


« 


Substituting  this  equation  in  Eq.  166,  the  optimality  criterion  can  be 
written  as 


m  _  Q.  . 

1  =  E  A-  —4—  (170) 

where 

«1J  -  *,(„>)  (171) 

and 

A.  =  - r—! -  (172) 

J  u  j 


Comparing  Eqs.  170  with  Eq.  65  shows  that  the  optimality  criterion  for 
the  system-stabi 1 i ty  constraint  problem  has  the  same  form  as  the 
optimality  criterion  for  the  displacement-constraint  problem.  Eq.  170 
can  also  be  written  as 


m 

1  =  l 

j=l 


q.  . 
A.  -^L 


J  p 


(173) 


l 


where 


’ij 


4i 


(174) 


is  the  energy  density  in  an  element  in  the  buckled  mode.  The  exponential 
recurrence  relation  and  the  linear  recurrence  relation  for  the  system- 
stability  constraint  problem  are  given  by 


and 


t(z 

’Vj-i  VX 


l/r 


(175) 


LL 


Z  A  .  - 

j-l  JP(e1A( 


-  1 


(176) 


Equations  to  evaluate  the  Lagrange  multipliers  a.  can  be  written  by  using 

•J 

the  procedure  discussed  for  the  multiple-displacement-constrained  problem. 
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SECTION  VII 


EFFECT  OF  STRUCTURAL  SYMMETRY  ON  THE  ALGORITHM 

In  the  previous  sections  when  deriving  the  optimality  criterion  and 
developing  the  algorithm,  we  had  not  taken  into  consideration  the  possible 
symmetry  of  the  structure.  The  symmetry  may  be  due  to  the  nature  of  the 
constraints  imposed  on  the  structure  and  also  due  to  the  multiple-loading 
conditions.  We  had  mentioned  that  the  multiple-loading  conditions  can  be 
treated  as  multiple  constraints.  A  substantial  amount  of  computational 
effort  can  be  saved  by  making  use  of  the  symmetry  of  the  structure  in  the 
analysis  phase  and  particularly  in  the  optimization  phase  of  the  algorithm. 
A  modification  to  the  optimization  algorithm  becomes  particularly  essential 
when  the  structure  is  idealized  with  a  large  number  of  elements,  the 
number  of  loading  conditions  is  more  than  one,  and  the  number  of  active 
constraints  is  large. 

We  will  consider  a  very  simple  case  and  illustrate  what  modifications 
to  the  algorithm  are  required.  Consider  a  symmetric  structure  subjected 
to  multiple-displacement  constraints.  Let  us  assume  that  there  are  five 
constraints  which  are  active  at  the  kth  iteration.  The  exponential  recur¬ 
rence  relation  can  be  written  as 


lk(  I  X.  Al 


077) 


Now  because  of  the  symmetry  of  the  constraints  let  us  assume  that 
X1  =x5=a1  anc*  *2=X4=X2  anc*  A3=A3*  "*^e  recurrence  relation  (Eq.  177)  can 


now  be  written  as 


au1  -  ak  f;  *  Qi5V  -  ( Qiz  *  M , ;  Y"  078) 


The  linear  equations  (Eq.  85)  to  determine  the  Lagrange  multipliers 
can  be  written  as 


^5x5  U}5  =  {C  }5 


(179) 
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where  Qt.  are  the  coefficients  of  the  square  matrix  multiplying  the  A 

vector  in  Eq.  85  and  C*  is  the  right  side  of  Eq.  85.  These  quantities 

are  evaluated  for  five  constraints.  In  Eq.  179  Of.  =  Q*. .  Replacing  the 

ij  ji 

A's  by  the  a's,  and  remembering  that  in  Eq.  179  because  of  symmetry  the 
first  equation  is  the  same  as  the  fifth  equation  and  the  second  equation 
is  the  same  as  the  fourth  equation,  we  can  write 

“l  ^11  +  ^15^  +  “2^12  +  ^1 4 ^  +  a3^13^  =  C1 

al  ^21  ^25^  *  a2^22  ^24 ^  +  a3^23^  =  ^2  (ISO) 

al  ^31  +  ^35^  +  a2^32  +  ^34 ^  +  a3^33^  =  C3 

These  three  equations  are  not  symmetric.  Taking  into  consideration  the 

★ 

symmetry  of  the  constraints  on  the  coefficients  Q..,  Eq.  180  can  be 

^  J 

written  as 

(Ql!  +  Q]5)ai  +  2Qi2a2  +  ^1 3a3  =  C1 

2Q]  2al  +  ^22  +  ^24  ^a2  +  Q23a3  =  C1  ^181^ 

2Qi3«i  +  2Q23a2  +  Q33a3  =  C3 

Thus  because  of  the  symmetry  we  have  to  solve  three  equations  instead  of 

* 

five.  In  addition  Eq.  179  would  require  evaluation  of  fifteen  Q. .  coef- 

'  J 

ficients  as  compared  with  eight  coefficients  for  Eq.  181.  For  this  small 
problem  this  modification  does  not  appear  to  be  a  substantial  reduction 
in  the  computational  effort.  However,  for  a  structure  with  a  large  num¬ 
ber  of  design  variables  and  potentially  active  constraints,  it  would  save 
considerable  CP  time.  The  modification  to  the  algorithm  discussed  here 
can  be  extended  to  other  types  of  constraints,  to  other  symmetric  condi¬ 
tions,  and  also  to  other  algorithms  discussed  previously. 

Note  that  if  in  a  structure  all  the  ' p '  Lagrange  multipliers 
associated  with  all  the  active  constraints  are  equal,  because  of  the  sym¬ 
metry,  all  ' p ’  equations  will  reduce  to  a  single  equation.  This  will 
require  evaluation  of  only  one  Lagrange  multiplier.  Thus  an  apparently 
multiple-constraint  problem  is  reduced  to  an  equivalent  single-constraint 
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problem.  For  such  a  problem  an  algorithm  based  on  the  single  dominant 
constraint  would  give  the  same  design  as  obtained  by  using  the  multiple- 
constraint  algorithm.  The  72-bar  tower  problem  considered  in  References 
1,  2,  29,  and  other  references  fall  under  this  category. 
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SECTION  VIII 
OPTIMIZATION  PROCEDURE 

In  Sections  III  through  VII  we  have  derived  the  optimality  criterion 
the  different  recurrence  relations  to  modify  the  design  variables,  and 
the  equations  to  estimate  the  Lagrange  multipliers.  We  will  discuss  in 
this  section  how  different  algorithms  can  be  developed  based  on  these 
equations  and  what  options  we  have  in  selecting  an  algorithm. 

The  main  steps  in  the  optimization  of  a  structure  are: 

(1)  Assign  initial  sizes  to  all  the  elements.  For  all  example 
problems  we  will  assume  them  equal. 

(2)  Analyze  the  structure  using  the  displacement  method  of  analysis 

(3)  Determine  the  stresses  in  the  elements. 

(4)  Scale  the  design  to  satisfy  all  the  constraints.  The  scaled 
design  is  a  feasible  design.  After  scaling,  at  least  one  constraint 
will  be  on  the  constraint  surface. 

(5)  Separate  the  constraints  into  potentially  active  and  passive 
categories. 

(6)  Separate  the  elements  into  active  and  passive  categories. 

(7)  Determine  the  flexibility  coefficients  corresponding  to  the 
potentially  active  constraints. 

(8)  Evaluate  the  Lagrange  multipliers. 

(9)  Modify  the  design  variables  using  a  recurrence  relation. 

(10)  Terminate  the  iterations  if  the  specified  criterion  is  satis¬ 
fied  or  go  to  step  2.  The  terminating  criterion  may  be  the  number  of 
iterations  or  the  percentage  difference  in  the  weight  of  the  structure 
between  the  consecutive  iterations. 
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The  first  four  steps  are  common  to  all  the  algorithms.  The  sub¬ 
sequent  steps  differ  depending  upon  (a)  the  division  of  active  and  passive 
constraints,  (b)  the  value  of  the  parameters  selected  to  control  the 
step  size,  (c)  the  method  used  to  determine  the  Lagrange  multipliers, 

(d)  the  recurrence  relation  used  to  modify  the  design  variables,  and  (e) 
the  number  of  times  the  design  variables  are  modified  before  reanalyzing 
the  structure.  We  will  discuss  in  detail  the  different  options  available 
under  these  categories. 

(a)  Active  and  Passive  Constraints 

The  active  constraints  or  potentially  active  constraints  are 
those  that  are  included  in  the  constraint  equations.  The  flexibility 
coefficients  and  the  Lagrange  multipliers  associated  with  the  active 
constraints  enter  the  optimality  criterion  and  the  recurrence  relation. 

The  remaining  constraints  which  are  imposed  on  the  structure  are  called 
passive  constraints,  and  these  are  satisfied  by  the  simple  approach  of 
scaling  the  individual  design  variables.  In  the  case  of  a  displacement 
and  stress  constraint  problem,  we  can  use  one  of  the  following  approaches 
to  design  a  structure. 

1)  Treat  all  the  displacements  as  well  as  the  stress  constraints 
as  potentially  active  constraints,  determine  the  Lagrange  multipliers 
associated  with  all  the  active  constraints,  and  use  them  in  the  recurrence 
relation.  In  this  case  only  the  elements  whose  sizes  are  governed  by 
minimum  or  maximum  size  constraints  would  be  passive. 

2)  Consider  only  the  displacement  constraints  as  potentially 
active  and  all  the  stress  constraints  as  passive.  Here,  the  stress  con¬ 
straints  are  not  directly  included  in  the  constraint  equations.  Elements 
in  which  the  stresses  are  greater  than  the  maximum  allowable  stress  are 
sized  by  using  the  FSD  algorithm.  Thus  the  size  of  an  element  is  a  maximum 
of  the  sizes  obtained  by  (1)  the  recurrence  relation  for  the  displacement 
constraints,  (2)  the  FSD  algorithm,  or  (3)  the  minimum  size  constraint. 

Any  element  whose  size  is  not  equal  to  the  one  obtained  by  using  a  recur¬ 
rence  relation  is  characterized  as  a  passive  element.  Since  the  number 

of  elements  and  which  elements  are  passive,  affect  the  equations  used  to 
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determine  the  Lagrange  multipliers  and  consequently  the  recurrence 
relation,  one  must  keep  track  of  the  passive  elements  and  update  the 
list  whenever  any  change  occurs. 

3)  Treat  only  one  or  some  dominant  displacement  constraints 

as  potentially  active  and  consider  the  remaining  displacement  constraints 
and  all  the  stress  constraints  as  passive. 

4)  Treat  all  the  displacement  constraints  as  being  independent 
of  each  other,  and  the  stress  constraints  as  passive  constraints.  In 
this  case  the  size  of  an  element  is  a  maximum  of  the  sizes  obtained  by 

(1)  the  recurrence  relation  for  the  single  displacement  constraint  applied 
to  each  displacement  constraint  separately,  (2)  the  FSD  algorithm,  or 
(3)  the  minimum  size  constraint.  This  method  is  known  as  the  envelope 
method. 

For  a  structure  subjected  to  stress  constraints  only,  we  can  use 
one  of  the  following  approaches: 

1)  Treat  all  the  stress  constraints  as  potentially  active, 
determine  the  Lagrangian  multipliers  associated  with  the  active  constraints, 
and  use  them  in  the  recurrence  relation.  The  passive  elements  in  this 

case  will  be  the  elements  whose  sizes  are  governed  by  the  minimum  size 
limits. 

2)  Treat  all  the  stress  constraints  as  active,  but  approximate 
the  virtual  load  system  with  the  actual  load  system.  In  this  case  the 
Lagrange  multipliers  are  not  determined.  The  elements  are  sized  on  the 
basis  of  the  strain  energy  stored  in  the  element  due  to  the  applied  load. 

3)  Treat  all  the  stress  constraints  as  active,  but  use  the 
FSD  algorithm. 

(b)  Step  Size 

The  convergence  behavior  depends  on  the  parameter  ' r '  in  the 
recurrence  relation  and  also  the  equations  used  to  determine  the  Lagrange 
multipliers  and  the  parameter  'b'  in  the  Newton-Raphson  iterative  method. 

We  have  derived  a  relationship  between  b,  <f>,  and  r.  For  structural 
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optimization,  r= 2  is  generally  adequate.  However  for  some  problems, 

depending  on  the  behavior  of  the  constraints,  'r'  must  be  increased  to 

reduce  the  step  size  and  get  better  convergence.  The  parameter  'r‘  can 

be  kept  constant  through  all  the  iterations  or,  based  on  certain  criteria, 

can  be  changed  after  each  iteration.  A  good  indication  of  the  proper 

selection  of  ' r '  is  a  reduction  in  the  weight  of  the  structure  after 

each  iteration.  Another  approach  in  controlling  the  step  size  is  to 

generate  an  intermediate  design  vector  whenever  the  weight  of  structure 

is  more  than  the  previous  iteration.  This  can  be  done  by  taking  the 

average  of  the  design  variables  in  the  previous  and  the  present  iteration. 

For  example,  if  the  weight  of  the  structure  with  the  design  variables 
k+1  k 

A*  is  greater  than  the  one  with  Ai  ,  the  intermediate  design  variable 
would  be  (a|+^  +  A^)/2.  This  process  can  be  continued  until  one  obtains 
a  design  with  a  weight  less  than  the  weight  of  the  previous  iteration. 

(c)  Lagrange  Multipliers 

When  the  structure  is  subjected  to  inequality  constraints,  we 
use  only  those  Lagrange  multipliers  which  are  positive.  A  negative  Lagrange 
multiplier  is  permissible  only  if  the  specific  constraint  is  to  be  enforced 
as  an  equality  constraint.  The  Lagrange  multipliers  can  be  determined 
by  one  of  the  following  methods: 

1)  Linear  simultaneous  equation. 

Eqs.  86,  87  for  the  displacement  constraints. 

Eq.  136  for  the  stress  constraints. 

2)  Newton-Raphson  iteration. 

Eq.  101  for  the  displacement  constraints. 

Eq.  141  for  the  stress  constraints. 

3)  Exponential  aid  linear  recurrence  relations. 

Eqs.  74  and  76  for  the  displacement  constraints. 

Eqs.  132  and  133  for  the  stress  constraints. 

4)  Approximate  relations. 

Eq.  95  for  the  displacement  constraint. 

Eqs.  138,  139,  140  for  the  stress  constraints. 
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Under  the  last  category  we  can  include  a  number  of  degenerate  equations 
obtained  from  the  first  method.  For  the  first  two  cases  and  some  equations 
of  the  fourth  case,  the  coefficients  must  be  calculated  Q, .  or  RiS  before 
the  Lagrange  multipliers  can  be  determined.  In  the  third  case  the  Lagrange 
multipliers  are  estimated  based  on  the  ratio  of  the  actual  value  of  the 
constraint  to  the  limiting  value  of  the  constraint,  and  the  coefficients 
or  F^j  need  not  be  known  before  the  x's  are  estimated.  The  advantages 
and  the  disadvantages  of  using  these  methods  were  discussed  in  the  section 
where  they  were  derived. 

(d)  Recurrence  Relations 

In  the  case  of  the  displacement  constraint  problem  there  are 
basically  three  distinct  relations.  These  are: 

1)  Exponential  relation  for  A^  or  zi  (Eq.  70  or  156). 

2)  Linear  relation  for  A^  (Eq.  71  or  72). 

3)  Linear  relation  for  z1-  (Eq.  157  or  158). 

.  A  tu 

We  can  use  any  one  of  these  relations.  In  these  relations  we 
can  keep  the  step-size  parameter  ' r '  constant  or  change  it  after  each 
iteration.  All  these  relations  need  the  evaluation  of  the  Lagrange 
multipliers. 

For  the  stress  constraint  problem  we  can  use  one  of  the  following: 

1)  Exponential  relation  for  A^  or  z^  (Eq.  130  or  162). 

2)  Linear  relation  for  A.  (Eq.  131). 

3)  Linear  relation  for  z^  (Eq.  163). 

4)  FSD  algorithm  (Eq.  117). 

5)  Exponential  relation  for  Ai  based  on  the  strain  energy 
distribution  (Eq.  147). 

6)  Linear  relation  for  A^  based  on  the  strain  energy  dis¬ 
tribution  (Eq.  148). 

The  first  three  of  these  relations  require  evaluation  of  the  Lagrange 
multipliers  before  they  can  be  used.  The  fourth  requires  knowing  the 
force  in  the  element,  and  the  last  two  use  the  strain  energy  stored  in 
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an  element.  For  the  last  three  the  Lagrange  multipliers  are  not  evaluated, 
only  the  information  from  the  applied  load  vector  is  used. 

e)  The  flexibility  coefficients  Q^-  used  to  evaluate  the  Lagrange 
multipliers  and  modify  the  design  variables  are  calculated  with  the  assumed 
values  of  the  design  variables.  When  we  modify  the  design  variables, 
the  previously  calculated  flexibility  coefficients  are  no  longer  valid. 

We  have  two  choices.  One  is  to  reanalyze  the  structure  and  calculate 
the  new  flexibility  coefficients.  And  the  other  choice  is  to  assume 
that  the  flexibility  coefficients  are  not  substantially  affected  with 
small  changes  in  the  design  variables  and  reuse  them.  The  second  option 
would  require  a  smaller  number  of  analyses.  However,  note  that  the  design 
does  not  move  away  from  the  region  where  the  flexibility  coefficients 
are  no  longer  valid. 
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SECTION  IX 

ILLUSTRATIVE  PROBLEMS 

This  section  presents  results  for  the  sample  problems  which  are 
solved  by  using  some  of  the  algorithms  and  the  optimization  procedure 
discussed  in  the  previous  sections.  Some  of  the  problems  selected  for 
presentation  are  standard  problems  and  have  been  used  by  several  inves¬ 
tigators  to  study  the  performance  of  their  algorithms.  The  first  problem 
is  a  10-bar  truss.  This  structure  is  designed  to  satisfy  stress- 
displacement  constraints  and  different  allowable  stresses  in  the  elements. 
The  second  structure  is  a  200-bar  truss  subjected  to  5  loading  conditions. 
This  structure  is  designed  with  constraints  on  stresses  and  displacements. 
The  third  structure  is  a  cantilever  box-beam  with  the  top  and  bottom 
skins  made  of  a  layered  composite  material.  The  box-beam  is  designed  to 
satisfy  stress  constraints  and  local  plate  buckling.  The  fourth  structure 
is  a  composite  wing  structure  designed  to  satisfy  stress  constraints  in 
all  the  elements  and  a  twist  constraint  at  the  free  end.  The  last  structure 
is  a  three-dimensional  dome  designed  to  satisfy  the  constraint  on  system 
stability.  For  all  the  problems,  unless  otherwise  mentioned,  the  sizes 
of  all  the  elements  are  equal  for  the  first  iteration.  The  CP  time  given 
for  some  problems  is  for  the  CDC  CYBER  74/175. 

In  problems  where  the  linear  equations  (Eqs.  86  or  136)  are  solved 
to  determine  the  Lagrange  multipliers,  the  equations  were  first  arranged 
in  the  order  of  the  degree  of  activity  of  the  constraints.  The  first 
equation  corresponds  to  the  constraint  that  has  been  closest  to  the  con¬ 
straint  surface,  and  the  last  equation  corresponds  to  the  constraint 
that  was  farthest  from  the  constraint  surface.  The  rearranged  equations 
were  solved  by  an  iterative  method.  In  the  first  iteration,  only  one 
equation  was  solved,  assuming  that  there  was  only  one  constraint,  and 
during  the  subsequent  iterations  a  new  equation  was  added  each  time. 

This  process  was  continued  until  all  the  equations  corresponding  to  all 
the  potentially  active  constraints  were  solved. 
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EXAMPLE  1.  TEN-BAR  TRUSS  (DISPLACEMENT  AND  STRESS  CONSTRAINTS) 

The  truss  shown  in  Figure  1  was  designed  to  satisfy  a  stress  limit 
of  25  ksi  in  all  the  elements  and  a  displacement  limit  of  +2.0  in.  at 
all  the  node  points  in  the  x  and  y  direction. 

Case  I 

The  truss  is  designed  by  considering  stress  constraints  as  well  as 
displacement  constraints  as  potentially  active  constraints.  The  three 
recurrence  relations  used  to  modify  the  design  variables  were  (1)  Case 
la,  the  exponential  relation  for  A^  or  zi  (Eqs.  70,  130,  156,  162),  (2) 
Case  lb,  the  linear  relation  for  A.  (Eqs.  71  and  131),  (3)  Case  Ic,  the 
linear  relation  for  z,-  (Eqs.  157  and  163).  The  Lagrange  multipliers  for 


E»10  psi 


p=0.1  lb/in” 


Min.  Size 


Figure  1.  Ten-Bar  Truss 
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the  three  cases  are  determined  by  solving  the  linear  simultaneous  equations 
given  in  Eqs.  86  and  136.  The  potentially  active  constraints  were  deter¬ 
mined  at  each  iteration.  For  the  first  iteration,  all  the  constraints 
satisfying  the  conditions  Cj/“Cj  =  s  =0.5  for  the  scaled  design  were 
assumed  to  be  potentially  active.  The  maximum  value  of  is  equal 

to  unity  for  the  scaled  design.  In  subsequent  iterations,  Cj/Tj  was 
determined  based  upon  the  activity  of  the  constraints  in  the  previous 
iteration.  For  example,  6  for  the  (k+l)th  iteration  is  equal  to  Cj/Cj 
of  the  kth  iteration  that  includes  all  the  active  constraints.  In  addition 
a  check  was  made  to  make  sure  that  there  was  at  least  one  nonactive  con¬ 
straint  included  in  the  potentially  active  constraints.  This  was  essential 
in  order  to  assure  that  all  the  active  constraints  were  used  in  the  recur¬ 
rence  relation.  The  step-size  parameter  'r'  for  all  three  cases  was  set 
equal  to  2.  The  structure  was  analyzed  after  each  iteration. 

The  iteration  history  for  Cases  la  through  Ic  is  given  in  Table  1. 

The  table  also  contains  the  active  constraints  at  each  iteration  and  the 
total  time  needed  to  complete  the  iterations  for  Case  Ic.  For  Case  la, 
after  the  12th  iteration,  the  weight  of  the  structure  was  found  to  oscil¬ 
late  at  weights  greater  than  5070.1  lb.  This  behavior  was  due  to  the 
large  step  size.  If  the  step  size  parameter  'r'  had  been  increased  after 
the  12th  iteration,  in  order  to  reduce  the  step  size,  the  additional 
iterations  would  lead  to  a  design  with  a  weight  of  5060.85  lb.  The  active 
constraints  at  the  optimum  for  this  design  were  the  vertical  displacement 
at  node  1  and  the  stress  in  member  5.  For  Case  la,  at  the  10th  iteration, 
the  scaled  weight  of  the  structure  jumped  to  6069.3  lb,  because  the  stress 
In  element  5  suddenly  became  active.  The  iterations  for  Cases  lb  and  Ic 
lead  to  a  design  with  a  minimum  weight  of  5076.6  lb  with  the  vertical 
displacements  at  nodes  1  and  3  being  active  and  equal  to  2.0  in.  The 
stress  constraint  for  Cases  lb  and  Ic  did  not  become  active.  Both  of 
the  designs  satisfy  the  optimality  criterion.  The  design  with  a  weight 
of  5076.6  lb  is  a  relative  minimum.  The  iteration  history  for  the  three 
cases  is  shown  in  Figure  2.  The  computer  time  shown  in  this  figure  is 
for  Case  Ic.  The  total  time  required  to  complete  12  Iterations  for  Case 
la  and  14  iterations  for  Case  lb  was  1.78  and  1.88  seconds,  respectively. 
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Figure  2.  Iteration  History  for  10-Bar-Truss 
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Case  II 

The  truss  is  designed  by  using  the  expdhential  recurrence  relation 
with  the  step-size  parameter  ' r *  set  equal  to  4  after  each  analysis  of 
the  structure.  The  Lagrange  multipliers  were  determined  by  using  Eqs. 

86  and  136.  For  Case  1 1  a,  the  structure  was  reanalyzed  after  each 
iteration.  For  Case  I I b ,  a  maximum  of  10  subiterations  were  allowed 
before  reanalyzing  the  structure  with  a  criterion  set  to  satisfy  the 
constraints  equal  to  10“^.  During  the  subiterations  the  step-size 
parameter  was  doubled,  to  reduce  the  step  size,  whenever  the  scaled 
weight  of  the  structure  was  greater  than  the  previous  lower-weight 
design  within  the  subiterations.  We  have  referred  to  this  approach  in 
the  last  seciton  as  the  modified  Newton-Raphson  approach.  The  iteration 
history  and  the  total  time  needed  to  complete  the  specified  number  of 
iterations  are  given  in  Table  2.  Case  1 1  a  and  Case  lib  gave  minimum 
weight  designs  of  5076.6  lb  and  5060.8  lb  respectively.  However,  the 
approach  of  using  repeated  subiterations  (Case  lib)  needed  substantially 
more  computer  time. 


TABLE  2 

ITERATION  HISTORY  OF  10-BAR-TRUSS 
STRESS-DISPLACEMENT  CONSTRAINTS  -  CASE  II 


Case  I  la 

Case  I lb 

Iteration 

No. 

Weight 

Total  Time 

In  Seconds 

Weight 

Total  Time 

In  Seconds 

1 

8266.1 

0.09 

8266.1 

0.43 

2 

6423.2 

0.21 

6017.7 

0.92 

3 

6009.5 

0.34 

5818.9 

1.39 

4 

5856.1 

0.45 

5689.0 

1.82 

5 

5780.1 

0.58 

5572.6 

2.36 

6 

5719.4 

0.70 

5445.3 

2.90 

7 

5660.4 

0.81 

5303.1 

3.50 

8 

5599.3 

0.93 

5198.9 

4.16 

9 

5536.6 

1.04 

5105.3 

4.95 

10 

5469.9 

1.16 

5079.6 

5.59 

11 

5398.4 

1.23 

5074.3 

6.59 

12 

5320.4 

1.41 

5064.0 

7.21 

13 

5232.3 

1.56 

5062.4 

7.96 

14 

5155.0 

1.69 

5061.2 

8.17 

15 

5076.9 

1.84 

5061.0 

9.46 

16 

5076.6 

2.00 

5060.8 

10.20 

57 
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The  iteration  history  for  Case  II  is  shown  in  Figure  3.  The 
Lagrange  multipliers  and  the  areas  of  the  members  for  the  two  minimum- 
weight  designs  are  given  in  fable  3,  and  are  designated  as  Design  1  and 
Design  2. 

EXAMPLE  2.  TEN-BAR  TRUSS  (STRESS  CONSTRAINTS) 

The  truss  shown  in  Figure  1  was  designed  to  satisfy  different 
stress  constraints  in  the  elements. 

Case  I 

The  maximum  allowable  stress  for  all  the  elements  was  25  ksi 
except  the  maximum  allowable  stress  was  increased  to  50  ksi  for  element 
9.  The  FSD  design  for  this  case  is  not  the  correct  minimum-weight 
design.  The  minimum  weight  obtained  by  FSD  algorithm  was  1725  1b. 

The  Lagrange  multipliers  for  this  case  were  obtained  by  solving  Eq. 
136  and  all  10  stress  constraints  were  assumed  to  be  potentially  active 
in  all  the  iterations.  The  structure  was  reanalyzed  after  each 
iteration.  The  parameter  'r'  was  set  equal  to  2  for  all  three 
recurrence  relations:  the  exponential  relation  (Eq.  130),  the  linear 
relation  for  A^  (Eq.  131)  and  the  linear  relation  for  zn-  (Eq.  163).  The 
parameter  r=2  was  found  to  be  too  large  to  obtain  proper  convergence  for 
the  exponential  relation.  Therefore,  these  results  are  not  given  here. 
The  iteration  history  for  the  other  two  cases.  Case  la  (Eq.  131)  and 
Case  lb  (Eq.  163),  is  given  in  Table  4.  This  table  contains  the  active 
constraints  at  each  iteration  and  the  total  time  required  to  complete 
the  specified  number  of  iterations.  The  iteration  history  is  shown  in 
Figure  4.  The  minimum-weight  design  for  this  problem  is  1497.6  lb  with 
the  stresses  In  elements  1,  2,  3,  4,  6,  7,  8,  10  equal  to  25  ksi.  The 
stress  at  the  optimum  in  element  9  is  37.5  ksi,  and  this  is  not  an 
active  constraint.  The  cross-sectional  areas  of  all  the  elements  and 
the  Lagrange  multipliers  associated  with  the  minimum-weight  design  are 
given  in  Table  3  as  Design  3. 


Iteration  History  for  10-Bar-Truss 
Stress-Displacement  Constraints  -  Case  II 


MINIMUM-WEIGHT  DESIGNS  OF  10-BAR  TRUSS 


TABLE  4 
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Figure  4.  Iteration  History  for  10-Bar-Truss 
Stress  Constraints  -  Case  I 
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Case  II 

The  problem  solved  was  the  same  as  Case  I  above,  except  that  the 
Lagrange  multipliers  were  determined  by  using  the  recurrence  relation  in 
Eq.  132.  The  step-size  parameter  r  was  set  equal  to  2,  and  the  parameter 
' b '  in  Eq.  132  was  set  equal  to  0.5  (See  Eq.  98).  The  three  recurrence 
relations  used  to  modify  the  design  variables  were  the  exponential  relation 
(Eq.  130),  the  linear  relation  for  Ai  (Eq.  131)  and  the  linear  relation 
for  zi  (Eq.  163).  The  Lagrange  multipliers  for  the  first  iteration  were 
assumed  to  be  proprotional  to  the  forces  in  the  bars.  All  Lagrange  multi¬ 
pliers  were  normalized  so  that  the  maximum  value  of  the  Lagrange  multiplier 
was  equal  to  unity.  The  structure  was  reanalyzed  each  time  the  design 
variables  were  modified  using  the  recurrence  relations.  The  exponential 
relation  gave  a  minimum  weight  design  of  1500.6  lb  after  13  iterations. 

(See  Table  5.)  Use  of  both  linear  relations  caused  oscillations  in  the 
weight  of  the  structure  and  did  not  give  a  good  minimum-weight  design. 
Therefore  the  problem  was  resolved  by  using  the  three  recurrence  relations. 
However,  in  the  first  iteration  the  design  variables.  However,  in  the 
first  iteration  the  design  variables  were  modified  by  using  the  FSD  algorithm 
(Eq.  117).  The  iteration  history  for  the  three  cases,  1 1  a  (Eq.  130), 
lib  (Eq.  131),  and  lie  (Eq.  163)  is  given  in  Table  5.  The  table  also 
contains  the  total  computer  time  required  to  complete  the  specified  number 
of  iterations.  This  time  is  substantially  less  than  that  given  in  Table 
4,  where  Eq.  136  was  used  to  determine  the  Lagrange  multipliers.  The 
iteration  history  for  the  three  cases  is  shown  in  Figure  5.  Previous 
experience  with  Eq.  132  has  shown  that  with  the  exponential  recurrence 
relation  with  r=2  and  the  parameter  b  increased  to  4  to  reduce  the  step 
size,  a  design  with  a  weight  1497.6  lb  can  be  obtained.  However,  this 
needs  a  substantially  larger  number  of  iterations. 

Case  III 

The  10-bar  truss  with  stress  constraints  was  solved  with  a  maximum 
allowable  stress  in  element  9  equal  to  37.5  ksi  and  a  maximum  allowable 
stress  in  all  other  elements  equal  to  25  ksi.  The  stress  in  element  9 
was  selected  on  the  basis  that  its  stress  does  not  increase  beyond  37.5 
ksi  for  the  minimum-weight  design.  The  structure  was  designed  by  using 
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Figure  5.  Iteration  History  for  10-Bar-Truss 
Stress  Constraints  -  Case  II 
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two  approaches.  In  Case  Ilia,  the  FSD  algorithm  (Eq.  117)  was  used. 

And  in  Case  Illb,  the  design  variables  were  modified  by  using  the  exponential 
relation  (Eq.  130)  and  the  Lagrange  multipliers  were  determined  by  using 
Eq.  140.  This  equation  assumes  that  there  is  no  coupling  effect  and 
that  the  elements  are  fully  stressed.  The  iteration  history  for  the  two 
cases  and  the  computer  time  for  Case  Illb  are  given  in  Table  6.  The  FSD 
algorithm  gives  a  fully  stressed  design  which  is  nonoptimum.  The  use  of 
Eq.  130  to  determine  the  Lagrange  multipliers  for  Case  Illb  gave  a  design 
with  a  weight  of  1502.8  lb  which  is  fully  stressed  but  closer  to  the 
real  optimum  weight  of  1497.6  lb.  If  the  allowable  stress  in  element  9 
was  increased  to  50  ksi,  the  use  of  Eq.  130  gives  a  design  with  a  weight 
of  1621  lb,  which  is  in  between  the  optimum  design  (1497.6  lb)  and  the 
design  (1725  lb)  obtained  by  using  the  FSD  algorithm.  The  iteration 
history  for  Case  Ilia  and  Illb  is  also  given  in  Figure  6. 

TABLE  6 

ITERATION  HISTORY  FOR  10-BAR-TRUSS 
-  STRESS  CONSTRAINTS  -  CASE  III 


Case  Ilia 

Case  Illb 

Total  Time 
in  Seconds 

3434.9 

3434.9 

msm 

1815.5 

1787.9 

HUS' 

1726.2 

1667.6 

■  1 

1648.1 

1569.2 

1590.6 

1545.4 

0.40 

1570.4 

1539.2 

0.48 

1568.5 

1531.1 

0.55 

1523.9 

0.63 

1517.3 

0.70 

1512.6 

0.78 

1504.8 

0.86 

1502.0 

0.94 
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EXAMPLE  3.  TWO-HUNDRED-BAR  TRUSS  (DISPLACEMENT  AND  STRESS  CONSTRAINTS) 

The  200-bar  truss  shown  in  Figure  7  is  subjected  to  five  loading 
conditions  as  follows: 

(1)  A  load  of  1000  lb  in  the  positive  X  direction  applied  at  nodes 

1,  6,  15,  20, . ,  71. 

(2)  A  load  of  1000  lb  in  the  negative  X  direction  applied  at  nodes 

5,  14,  19, . ,  75. 

(3)  A  load  of  1000  lb  in  the  negative  Y  direction  at  nodes  1,  2, 

3,  4,  5,  6,  8,  10  12,  14,  15,  16,  . ,  73,  74,  75. 

(4)  Loading  conditions  1  and  3. 

(5)  Loading  conditions  2  and  3. 

The  maximum  allowable  stress  in  all  the  elements  is  equal  to  +10 
ksi.  The  displacement  limit  in  the  X  and  Y  directions  at  all  the  nodes 
is  +0.5  in. 

The  structure  was  designed  to  satisfy  stress  and  displacement  con¬ 
straints.  Two  approaches  were  used.  In  the  first  approach  the  constraints 
on  the  displacements  and  the  stresses  were  both  treated  as  potentially 
active  constraints,  i.e.,  the  Lagrange  multipliers  associated  with  the 
stress  as  well  as  the  displacement  constraints  were  determined.  In  the 
second  approach,  the  stress  constraints  were  treated  as  passive  constraints, 
i.e.  no  Lagrange  multipliers  or  gradients  of  the  constraints  associated 
with  stresses  in  the  elements  were  evaluated.  The  Lagrange  multipliers 
corresponding  to  only  the  displacement  constraints  were  used  in  the  recur¬ 
rence  relation.  The  linear  recurrence  relations  for  A^  (Eqs.  71  and 
131)  were  used.  The  effect  of  using  other  recurrence  relations  may  be 
found  in  Reference  31. 

The  step-size  parameter  ' r '  was  set  equal  to  '2',  and  the  Lagrange 
multipliers  were  determined  by  using  the  linear  equations  (Eqs.  86,  136). 

The  structure  was  reanalyzed  for  both  cases  after  each  iteration.  The 
symmetry  of  the  structure  was  not  taken  into  consideration  in  order  to 
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reduce  the  loading  conditions  or  the  number  of  design  variables.  The 
number  of  potentially  active  constraints  is  determined  by  using  the  method 
discussed  in  Example  1.  The  iteration  history  for  the  two  cases,  Case  I 
(stress  +  displacements)  and  Case  II  (displacements),  is  given  in  Table 
7.  The  table  also  gives  the  CP  time  required  for  each  iteration.  The 
average  CP  time  to  analyze  the  structure  was  found  to  be  2.4  seconds. 

The  difference  between  the  time  given  in  the  table  and  2.4  seconds  is 
the  time  required  to  complete  an  optimization  phase  of  the  iteration. 

The  optimization  phase  includes  evaluating  the  constraint  gradients, 
assembling  the  elements  of  Eqs.  86  or  136,  solving  the  equations  to  deter¬ 
mine  the  Lagrange  multipliers,  and  modifying  the  design  variables  by 
using  the  recurrence  relation.  The  iteration  history  for  the  two  cases 
is  also  given  in  Figure  8. 

The  weight  of  the  optimum  design  for  this  structure  is  28858.4  lb. 

The  areas  of  the  elements  associated  with  this  design  may  be  found  in 
Reference  31. 

EXAMPLE  4.  CANTILEVER  BOX-BEAM 

The  structure  shown  in  Figure  9  was  idealized  with  membrane  quadri¬ 
lateral  elements,  shear  panels,  and  bar  elements  (posts).  The  quad¬ 
rilateral  membrane  elements  in  the  top  and  bottom  skins  consisted  of 
four  layers  with  fibers  in  the  0°,  90°,  +45°,  and  -45°  direction.  The 
0°  fibers  are  parallel  to  the  length  of  the  beam.  The  idealized  struc¬ 
ture  consisted  of  18  quadrilaterals,  18  shear  panels  and  18  posts.  The 
six  quadrilaterals  near  the  tip  of  the  box  beam  consisted  of  four  layers 
of  graphite  epoxy.  The  six  quadrilateral  elements  in  the  middle  con¬ 
sisted  of  boron  epoxy  in  the  0°  direction  and  the  remaining  three  layers 
of  graphite  epoxy.  The  six  quadrilateral  elements  closer  to  the  fixed 
end  consisted  of  all  layers  of  boron-epoxy.  The  shear  panels  and  posts 
have  different  material  properties.  The  composite  elements  were  designed 
by  using  the  maximum  stress  criteria.  The  material  properties  used  for 
this  example  are  given  in  Table  8.  The  reason  for  selecting  diverse 
properties  for  the  different  elements  was  to  illustrate  the  versatility 
of  the  computer  program.  The  program  used  to  solve  the  problem  is  docu¬ 
mented  in  Reference  33. 


TABLE  7 
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The  box  was  subjected  to  two  loading  conditions.  In  the  first 
loading  condition  1000  lb  was  applied  at  nodes  1  through  4  in  the 
negative  z  direction.  In  the  second  loading  condition,  500  lb  was 
applied  at  nodes  2  and  4  in  the  negative  y  direction.  The  elements  were 
designed  by  using  the  recurrence  relation  (Eq.  147)  based  on  the  strain 
energy  stored  in  the  element.  To  take  into  consideration  the  effect  of 
the  two  loading  conditions,  the  sum  of  the  strain  energies  for  both 
loading  conditions  was  used.  In  the  first  iteration  the  number  of 
layers  in  all  four  fiber  directions  for  all  the  elements  was  assumed  to 
be  equal.  The  beam  was  first  designed  to  satisfy  the  stress 
constraints.  Then  the  design  was  modified  to  prevent  local  buckling  of 
the  elements.  This  was  achieved  by  using  the  simple  expression  for  the 
buckling  of  simply  supported  orthotropic  plates  (see  Reference  33).  The 
thickness  of  the  plate  elements  was  increased  to  satisfy  the  buckling 
equations,  and  the  whole  structure  was  reanalyzed  to  determine  the 
modified  in-plane  forces.  This  process  was  continued  until  there  was  no 
further  increase  in  the  thickness  of  the  plate  elements.  Since  the 
buckling  load  of  the  plate  can  be  further  increased  by  increasing  the 
thickness  of  the  +45°  layers  instead  of  0°  or  90°  layers,  only  the 
thickness  of  these  layers  was  primarily  increased.  The  convergence 
behavior  of  the  box-beam  is  given  in  Table  9.  The  minimum  weight  of  a 
design  satisfying  only  the  stress  constraints  was  26.45  lb,  and  the 
design  with  stress  and  local  buckling  constraints  was  34.40  lb.  The 
distribution  of  the  number  of  plies  in  the  0°,  90°,  and  +45°  direction 
is  given  in  Figures  10  and  11.  The  number  of  plies  in  the  +45°  and  -45° 
direction  is  equal.  The  number  of  plies  in  each  layer  was  obtained  by 
dividing  the  thickness  of  each  layer  by  the  ply  thickness  and  rounding 
the  result  to  the  next  higher  integer.  The  number  of  plies  in  the  top 
skin  for  both  cases  is  nearly  the  same,  but  the  number  of  plies  in  the 
bottom  skin,  which  is  subjected  to  compressive  loads,  has  substantially 
changed.  The  change  in  the  number  of  plies  in  the  0°  and  90°  direction 
is  due  to  scaling  the  design  variables  and  rounding  the  thickness  to 
make  it  equal  to  a  multiple  of  a  single-ply  thickness. 
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TABLE  9 

ITERATION  HISTORY  FOR  CANTILEVER  BOX-BEAM 


Iteration 

No. 

Weight 

Iteration 

No. 

Weight 

1 

153.20 

10 

29.27 

2 

74.83 

11 

26.97 

3 

32.71 

12 

26.50 

4 

30.48 

13 

26.46 

5  ! 

29.27 

14 

26.45 

6 

29.34 

15 

32.08* 

7 

29.62 

16 

33.90* 

8 

29.56 

17 

34.26* 

9 

29.50 

18 

34.36* 

19 

34.40* 

*  -  modified  for  local  buckling 

C.P.  time  for  each  Iteration  (1  -  14)  3.5  seconds 

C.P.  time  for  each  iteration  (15  -  18)  4.9  seconds 


In  the  optimum  design  for  the  stress-constraint  problem,  note  that 
the  distribution  of  the  number  of  plies  in  the  top  and  bottom  skin  near 
the  root  is  not  the  same.  The  distribution  of  layer  thicknesses  for  the 
first  10  iterations,  until  the  weight  reached  29.27  lb,  was  symmetric; 
but  with  additional  iterations  the  weight  was  reduced  to  26.45  lb  and 
the  symmetry  was  lost.  This  may  be  due  to  unsymmetric  loading  con¬ 
ditions. 

EXAMPLE  5.  WING  STRUCTURE 

The  finite  element  model  of  the  wing  structure  (Reference  37)  shown 
in  Figure  12  was  designed  to  satisfy  (1)  stress  constraints,  (2)  negative 
twist  constraints  (a  wash-out  condition)  and  (3)  positive  twist  constraints 
(a  wash-in  condition).  The  structure  was  subjected  to  a  single  loading 
condition.  The  coordinates  of  the  node  points  and  the  loads  applied  at 
each  node  are  given  in  Table  10.  The  structural  model  has  88  nodes  and 
158  members.  The  top  and  bottom  skins  were  idealized  by  membrane  quadri¬ 
laterals  and  triangles.  The  aluminum  spars  and  ribs  were  idealized  with 
shear  panels  and  posts.  The  skin  elements  consist  of  four  layers  with 
fibers  in  0°,  90°,  +45°,  and  -45°  directions.  The  0°  fibers  are 
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Figure  10.  Optimum  Design  (Stress  Constraints) 
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Figure  11.  Optimum  Design  (Stress  and  Local  Buckling  Constraints) 


Figure  12.  Model  of  Wing 
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TABLE  10 


COORDINATES  OF  NODES  AND  APPLIED  LOADS 


JOINT 

X 

• 

i 

-< 

-Z 

■m 

■sea 

F09CE-Z  1 

1 

70. 533 

90.000 

1.213 

o .  c  o  :• 

c .  c  o : 

29. 0T  0 

2 

70,633 

90.000 

-1.31J 

o.ccc 

o.ooc 

29. C  00 

3 

78. 167 

90.000 

1.500 

-25.3,300 

-  6  9  6  3  .  C  >3  C 

1133. Cr0 

<* 

78.167 

90.000 

-1.500 

2s: o.ooc 

KQ  WW 

1 1  3  3  .  C  C  0 

5 

85.5  0  3 

9  0 . 0  0  C 

1.213 

0.000 

jjgilMF  Wm 

90.933 

6 

85.500 

90.000 

-1.313 

0.000 

■k  hw 

50.930 

7 

92.633 

90.000 

1.125 

-9370.000 

nua«  fill 

11 30. CC  0 

e 

92. 633 

90.000 

-1.125 

9  8  7  3  .  0  0  0 

97 6  3.  3  CC 

11 3 3.  CC  0 

9 

63.  50  0 

90.000 

1.125 

2  L  5 . 0  0  C 

-7  2  6  3.33 C 

9  2  6 . 0  C  0 

1j 

63.500 

9  0  .  a  *  c 

-1.125 

-235.  -J  0  0 

7363. COu  • 

926. OOC 

11 

69. 686 

87.471 

1.349 

0.000 

0.000 

17»>.  OCO 

12 

69.636 

97.471 

-1.249 

0.0^0 

O.OOC 

175.003 

13 

76.097 

84.851 

1.586 

0.000 

0.000 

214. COO 

1-4 

76. 0  97 

84.351 

-l.see 

3. TOC 

t  '  r  f* 

c  •  «.  J  - 

214. roo 

15 

82.  7 4.6 

82.133 

1.427 

0  .  C  0  c 

0.  coo 

253. COO 

16 

82.74.6 

82.132 

-1.427 

G  .  O'  0  0 

L'  I  ['  w  l 

253.CT  0 

17 

89.647 

79.  312 

1.255 

-565  3.0  0  0 

?  3  c.  u  •  6  j  r 

1  0  2  3  . COD 

ie 

89.  6-7 

79.312  | 

-1.256 

5663.000 

*2320.300 

1023.000 

19 

57. 266 

77.669  ! 

1.275 

2313.000 

-946. DOG 

723.C03 

2  C 

57. 266 

77.669 

-1.276 

-2310. 000 

946. o:  0 

723.  0r  0 

21 

63. 992 

74.920 

1.532 

0  •  o  0  0 

0.000 

314. OCO 

22 

63.992 

74.920 

-1.532 

3.000 

o.c  oc 

314.000 

23 

70.962 

72.071 

1.756 

0.000 

0.000 

326. COC 

24 

70.96? 

72.071 

-1.799 

0.000 

0.C0C 

326. OCO 

25 

78.191 

69.116 

1.6  17 

5 . :  r  c 

o.  ore 

33?. COC 

26 

76.191 

69.116 

-1.617 

0.000 

0.300 

33  6.  D  00 

27 

85.692 

66.050 

1.424 

-4370.  me 

3 6 60. COC  i 

9r2.000 

28 

85.692 

6 6.  0  5  l 

-1.424 

4370.000 

-1660.000 

9C2.000 

29 

51. : 32 

65.339 

1.433 

lKt.uCC 

-713. COC 

646. COC 

3 ' 

51.032 

65.339 

-1.433 

-1743,000 

713. 00C 

646.000 

31 

55. 297 

62.369 

1.715 

3.000 

0.303 

34  * • c :  0 

32 

59.297 

62.369 

-1.715 

3.000 

o.OOC 

340.000 

33 

65. 826 

59.291 

2.0  12 

3.000 

O.OOC 

3  5  2 . 0  : 0 

3*. 

65. 326 

59.291 

-2.C12 

O.OOC 

o.ooc 

352. COO 

35 

73.635 

56.100 

1.8C7 

O.OOC 

o.cco 

365. CC0 

.  36 

73.635 

56.100 

—  1.807 

3.030 

0.300 

365.000 

37 

81. 73? 

52.75  7 

1.59C 

-425 O.OOC 

1740,  GC  0 

974. Cf  C 

38 

81.733 

52.787 

-1.55C 

425 O.OOC 

-1740. COC 

974. 0CC 

39 

44. 799 

53. DC  8 

1.5  87 

1823. OCO 

-743.000 

694. CCO 

<*: 

44.799 

53.008 

-1.587 

-1823. COC 

743. 0  0  0 

694. CCC 

41 

52.603 

49. ei e 

1.89? 

O.OOC 

0.  C  00 

36  5. CC  0 

42 

1  52.603 

49.918 

-i.ese 

o.ooc 

0.  COC 

265. CCO 

43 

60.691 

46.512 

2.225 

O.OOC 

O.OOC 

378. OC  3 

44 

60.691 

46.512 

-2.225 

O.OOC 

0  .  C  0  0 

376.000 

45 

69. : 79 

43.063 

1.997 

3  .  C  3  L 

o .  c  c  r 

392.000 

46 

69. L79 

43.063 

-1.957 

3.000 

o.CCC 

392. OCO 

47 

77.764 

39.525 

1.756 

-4443.000 

1623. ./DC 

i  c  5  o .  c :  c 

48 

77.  764 

39.525 

-1.756 

-1623,000 

1C  5  0.0  0  0 

49 

38 • 565 

43.676 

1.742 

gumal 

-773. cc: 

742.000 
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TABLE  10  (CONTINUED) 


JOINT 


39.565 
46.938 
<♦6.5:9 
55. 555 
55.555 
64.52  3 
6<*.523 
.73.  630 
73.633 
32. 331 
32.331 
41.214 
41.214 
50.1-20 
50.420 
59.967 
59.  967 
69. e76 
69. 676 
25.166 
25. 166 
35.593 
35.683 
46. 181 
46.181 
56.964 


14.17 
12. 3C 
12.30 

10.40 

13.40 
6.46 

8.4  6 
6.50 

6 . 5  C 
0.0? 
O.Ou 
0.30 
o .  *.  r 


54. COO 

54. c:  o 

66 . u  3 3 

66.000 


*  -2 

■B 

FORCE-Y 

FORCE-Z 

-1. 

742 

-1890. CCO 

773.000 

742.000 

2. 

C  82 

o .  o r  e 

o .  c :  c 

393.020 

-  2  * 

C  6  2 

c.occ 

O.COC 

353. OCC 

2. 

4  39 

j  «  j  .  c 

o.  cor 

4  C  4  •  C  ~  0 

-2. 

43  9 

0 . 0  c  c 

o.ooc 

404.000 

2. 

197 

C  •  C  rJ  L 

o.ccc 

420.002 

-2. 

187 

o.ooc 

o.ooc 

420.000 

1. 

922 

-4640.000 

1902. OCC 

1123.000 

-1. 

522 

4643. OOC 

-19C3.CCC 

1123.000 

1. 

696 

2293.020 

-937. OOC 

e  e  3 .  c  c  o 

-1 . 

e96 

-2293.00: 

937.000 

883. COO 

2. 

265 

0.000 

o .  :■  2  o 

413.000 

-2. 

265 

0.200 

o.ccc 

413. CCC 

2. 

£51 

O.OOC 

O.OOC 

391.230 

-2. 

651 

O.COC 

c.oco 

391. CCO 

2. 

376 

0.000 
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parallel  to  the  direction  of  the  middle  spar  as  defined  by  connecting 
nodes  4  and  84  in  the  top  skin  or  3  and  83  in  the  bottom  skin.  The 
elastic  constants  and  the  allowable  strengths  for  graphite-epoxy  and 
aluminum  are  given  in  Table  11.  The  maximum  stress  criteria  was  used  to 
design  the  composite  elements.  For  the  first  iteration  for  all  three 
cases,  the  relative  sizes  of  the  bar  elements  and  the  thicknesses  of  the 
plate  element  were  1.0  to  0.1,  respectively.  The  percentage  of  plies  in 
the  four  fiber  directions  for  the  initial  design  were  assumed  to  be 
equal . 

(a)  Stress-Constraint  Design 

The  structure  was  designed  by  using  the  recurrence  relation 
(Eq.  147)  based  on  the  strain  energy  in  each  layer  and  in  each  element. 
The  initial  weight  of  the  structure,  with  an  equal  percentage  of  plies 
in  the  four  fiber  directions,  was  312.84  lb,  and  the  minimum-weight 
design  obtained  after  14  iterations  was  45.45  lb.  The  iteration  history 
is  given  in  Table  12.  The  weight  of  the  structure  after  7  iterations 
was  found  to  be  45.63  lb.  The  design  was  scaled  after  each  iteration  to 
obtain  a  feasible  design.  The  distribution  of  the  number  of  plies  in 
the  four  fiber  directions  is  shown  in  Figure  13.  The  tip  deflections  at 
nodes  7  and  9  are  given  in  Table  13.  The  tip  twists  through  an  angle  - 
5.458°.  This  is  the  angle  between  the  lines  joining  nodes  9  and  7 
before  and  after  loading. 

(b)  Twist-Constraint  Designs 

The  wing  structure  was  designed  to  satisfy  two  twist  con¬ 
straints  to  illustrate  the  use  of  the  algorithm.  The  two  angles  selected 
for  this  are  -7°  (wash-out)  and  +2°  (wash-in)  (see  Figure  14).  The 
structure  was  designed  by  using  Eq.  48  for  a  single-displacement  con¬ 
straint  where  the  virtual  load  consisted  of  a  unit  couple.  The  stress 
constraints  in  the  elements  were  treated  as  passive  constraints.  The 
iteration  history  for  the  wash-out  and  wash-in  cases  is  shown  in  Figures 
15  and  16,  respectively.  After  each  iteration  the  structure  was  scaled 
to  satisfy  the  stress  and  twist  constraints.  The  two  designs  in  Figures 
15  and  16  are  referred  to  as  the  'stress  design'  and  the  'twist  design'. 
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TABLE  11 

PROPERTIES  AND  ALLOWABLE  STRESSES  FOR  WING  STRUCTURE 


Properties 

Graphite  Epoxy 

Aluminum 

(psi) 

18.5  x  106 

10.5  x  106 

E22  (psi) 

1.6  x  106 

10.5  x  106 

v12 

0.25 

0.3 

Shear  Modulus 

0.65  x  106 

4.038  x  106 

Density  (lbs/in^) 

0.055 

0.1 

Lamina  Thickness  (in) 

0.0052 

Allowable  Stresses  (ksi) 

Fx  (tension) 

139.0 

45 

Fx  (compression) 

86.0 

45 

Fy  (tension) 

- 

45 

Fy  (compression) 

- 

45 

F  u 

46.8 

25.9 

xy 

TABLE  12 

ITERATION  HISTORY  FOR  WING  STRUCTURE  STRESS  CONSTRAINTS 


NUMBER  OF  LAYERS  WITH  FIBER  ORIENTATIONS  IN 
0°,  90°,  +45®,  -45°  DIRECTIONS,  RESPECTIVE!  Y 


Figure  15.  Iteration  for  Twist  Constraints  (Wash-out 


STRESS  DESIGN 


Iteration  for  Twist  Constra 
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In  the  'stress  design'  the  stresses  in  all  the  elements  are  less  than 
orequal  to  the  maximum  allowable  stress,  and  in  the  'twist  design'  the 
twist  of  the  tip  of  the  wing  is  equal  to  the  prescribed  value. 

For  the  wash-out  case  the  lowest  design  weight  of  71.86  lb  was 
obtained  in  the  11th  iteration.  The  weight  after  the  4th  iteration, 

72.22  lb,  did  not  improve  significantly  with  additional  iterations.  The 
displacements  at  nodes  7  and  9  for  the  least-weight  design  are  given  in 
Table  14,  and  the  distribution  of  the  plies  in  each  element  is  given  in 
Figure  17. 

In  the  wash-in  case  the  twist  at  the  tip  for  the  initial  design  was 
negative.  This  design  could  not  be  scaled  to  achieve  the  +2°  twist. 

From  the  2nd  to  the  8th  iteration  the  twist  of  the  stress  design 
increased  until  it  reached  +2°.  In  the  subsequent  iterations  for  all 
cases,  the  twist  designs  were  acceptable  since  they  satisfied  the  stress 
constraints  also.  Figure  16  shows  that  the  weight  of  the  structure 
oscillates  and  does  not  reach  a  stable  position  in  the  design  space. 

The  least-weight  design,  with  a  weight  of  105.54  lbs  was  obtained  at  the 
21st  iteration.  The  oscillations  may  be  due  to  the  large  value  of  the 
step-size  parameter,  rs2,  selected  for  this  problem. 

The  displacement  at  nodes  7  and  9  for  this  design  are  given  in 
Table  15.  Figure  18  shows  the  distribution  of  the  plies  in  each  element. 
The  figure  shows  a  heavy  concentration  of  0°  plies  only  on  one  side  of 
the  wing.  A  design  with  an  even  distribution  of  plies  would  have  resulted 
for  both  cases  if  the  0°  plies  were  placed  at  an  angle  instead  of  parallel 
to  the  line  joining  nodes  3  and  83. 

EXAMPLE  6.  DOME  STRUCTURE 

The  dome  structure  shown  in  Figure  19  was  optimized  to  satisfy  system 
stability.  The  structure  was  subjected  to  a  load  of  1000  lb  applied  in 
the  vertical  direction  at  each  node  point.  The  structure  has  61  nodes 
and  132  members.  The  dome  is  supported  at  all  the  node  points  on  the 
boundary  and  all  the  joints  are  hinged.  The  structure  was  idealized 
with  bar  elements  carrying  axial  load  only. 
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TABLE  13 

DEFLECTIONS  AT  TIP  FOR  THE  STRESS-CONSTRAINT  DESIGN 


Node  7  (trailing  edge) 
Direction  (in) 


1 

Y 

Z 

-0.06 

-.19 

13.43 

Node  9  (leading  edge) 
Direction  (in.) 

Twist  (deg.) 

MB 

Y 

Z 

-.04 

-.16 

10.62. 

-5.45 

TABLE  14 


DEFLECTIONS  AT  TIP  FOR  THE 
TWIST-CONSTRAINT  DESIGN 
(WASH-OUT  7°) 


Node  9  (leading  edge) 
Direction  (in) 

Twist  (deg.) 

M 

Y 

Z 

-.08 

-.10 

8.07 

-7.00 

Node  7 

(trailing  edge) 

Direction  (in) 

H 

Y 

Z 

-.09 

-.15 

11.68 

TABLE  15 

DEFLECTIONS  AT  TIP  FOR  THE 
TWIST- CONSTRAINT  DESIGN 
(WASH-IN  2°) 


Node  7 

[trailing  edge) 

Direction  (in) 

X 

Y 

Z 

1  Node  9 

(leading  edge)  L 

Direction  (in) 

HIR 

Y 

Z 

Twist  (deg.) 


FIBER  ORIENTATIONS  IN 
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E  = 

Stress  Limit  = 
Special  wt  = 
Minimum  Size  = 


lO^psi 

75 ,000psi 

0.1 

0.02  (Relative'' 
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The  structure  was  designed  by  using  the  recurrence  relation  (Eq. 
175)  which  uses  the  strain  energy  in  an  element  in  the  buckled  mode 
associated  with  the  lowest  buckling  load  of  the  structure.  The  stress 
constraints  in  the  elements  were  treated  as  passive  constraints.  The 
structure  was  designed  by  using  two  values  of  the  step-size  parameter, 
r=2  (Design  I)  and  r=10  (Design  II).  The  iteration  history  for  the  two 
cases  is  given  in  Tables  16  and  17,  respectively.  The  minimum-weight 
designs  are  given  in  Table  18. 

In  these  tables  'Design  B'  is  the  scaled  design  satisfying  the 
stress  constraints  in  the  elements.  The  buckling  load  of  'Design  E'  is 
equal  to  the  applied  load.  For  all  the  designs  the  system  buckling  is 
the  active  constraint.  The  stress  constraint,  because  of  the  high 
allowable,  does  not  become  active. 

The  step-size  parameter  r=2  was  too  large  to  obtain  good 
convergence.  Therefore  the  step-size  was  reduced  by  generating  the 
intermediate  design  vector  whenever  the  scaled  weight  of  the  structure 
after  each  iteration  was  greater  than  the  previous  lower-weight  design. 
(See  discussion  in  Section  VIII).  The  intermediate  designs  obtained  by 
this  procedure  are  indicated  by  '**'  in  Table  16.  In  the  case  of  the 
step-size  parameter  r=10,  the  minimum  weight  design  was  obtained  after 
three  iterations.  The  weights  of  the  optimum  designs  obtained  for  the 
two  step-size  parameters  are  nearly  equal.  However,  Table  18  shows  that 
the  two  designs  are  not  identical. 
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TABLE  16 

ITERATION  HISTORY  OF  DOME  STRUCTURE 
FOR  m  =  2  AND  INTERMEDIATE  DESIGN 
VECTORS  (DESIGN  I) 


Design  B 

Design  E 

Cycle 

Wt 

Buckling 

B.Load 

Wt 

No. 

(lbs) 

Load 

Wt. 

(lbs) 

1 

16.583 

0.1877 

0.0113 

87.836 

2 

91.767 

0.6073 

0.0066 

151.086 

3 

19.503 

0.2679 

0.0137 

72.774** 

4 

45.307 

0.4551 

0.0100 

99.537 

5 

25.144 

0.3084 

0.0122 

81.511 

6 

21.708 

0.3198 

0.0147 

67.859 

♦Load  in  kips  over  each  node  point 

♦♦Weight  corresponding  to  intermediate  design  vectors 

2 

Initial  Design  -  Area  of  all  members  *  0.2160  in 


TABLE  17 

ITERATION  HISTORY  OF  DOME  STRUCTURE 
FOR  m  =  10  (DESIGN  II) 


Design  B 

Design  E 

Cycle 

Wt 

Buckling* 

B.Load 

wt 

No. 

(lbs) 

Load 

WtJ 

(lbs) 

1 

16.583 

mm 

0.0113 

2 

21.157 

3 

24.512 

■■ 

1 

67.059 

♦Load  in  kips  over  each  node  point 

Initial  Design  -  Area  of  all  members  =  0.2160  in2 
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SECTION  X 

SUMMARY  AND  CONCLUSIONS 

The  minimum-weight  design  must  satisfy  the  optimality  criterion 
derived  for  the  specific  type  of  constraints  imposed  on  the  structure. 

m  Fi  i 

The  criterion  states  that  a  quantity  ip  =  z  — - —  associated  with 

all  the  elements  in  the  structure  is  equal  to  unity.  The  quantity  F. . 

*  J 

depends  on  the  nature  of  the  constraints.  In  the  case  of  displacement, 
stress  and  stability  constraints  F.^  is  equal  to  Q^,  R^  and 
respectively  and  is  defined  in  Eqs.  10,  125  and  171.  F^  is  a  function 
of  the  energy  stored  in  the  element,  not  the  nature  of  the  energy 
depends  upon  the  type  of  constraints. 

The  algorithms  based  on  an  optimality  criterion  are  iterative  and  use 
the  recurrence  relation  derived  from  the  optimality  criterion  to  modify 
the  design  variables.  There  are  basically  three  distinct  relations 
which  we  have  derived  from  the  optimality  criterion.  These  are  (1)  the 
exponential  relation,  (2)  the  linear  relation  for  A^ ,  and  (3)  the  linear 
relation  for  Many  other  recurrence  relations  can  be  written  from 
these  thre  by  selecting  different  values  of  the  step-size  parameter  r. 
Each  of  the  recurrence  relations  contains  the  quantity  <|».  Use  of  the 
recurrence  relation  moves  the  initial  design  to  the  optimum  design 
satisfying  the  optimality  criterion.  When  the  condition  ip=l  is 
satisfied,  the  recurrence  relation  does  not  change  the  design  variables. 
However,  In  the  optimization  algorithm  based  on  the  optimality  criterion 
some  of  the  constraints  are  treated  as  passive  and  the  step-size 
parameter  is  generally  kept  constant;  therefore  the  scaled  weight  of  the 
structure  reaches  a  certain  minimum  value,  and  then  it  increases.  The 
real  minimum  weight  lies  closer  to  the  lowest-weight  design  obtained  by 
using  the  algorithm.  The  real  minimum-weight  design  satisfying  the 
optimality  criterion  can  be  obtained  from  the  lowest  weight  with 
additional  iterations  but  with  a  smaller  step  size.  The  difference 
between  the  real  minimum  and  the  lowest-weight  design  is  generally 
small.  This  behavior  is  particularly  true  when  the  stress  dominate. 
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The  recurrence  relations  can  only  be  used  after  determining  the 
coefficients  and  the  Lagrange  multipliers.  The  coefficients  F^j  are 
related  to  the  gradients  of  the  constraints.  The  number  of  Lagrange 
multipliers  depends  upon  the  number  of  active  constraints.  The  active 
constraints  are  defined  as  those  that  are  associated  with  the  positive 
Lagrange  multipliers.  .  These  constraints  are  very  close  to  the  constraint 
surface.  The  number  of  active  constraints  changes  with  the  iterations. 
Initially  the  number  is  small,  but  as  the  design  approaches  the  optimum 
the  number  increases.  The  constraints  sometimes  switch  from  the  active 
category  to  the  passive  category,  depending  upon  their  activity  as  the 
optimum  process  proceeds  from  one  iteration  to  the  next.  A  prior  knowl¬ 
edge  of  the  active  constraints  at  the  optimum  is  not  necessary.  At  each 
iteration  one  has  to  determine  the  set  of  active  constraints. 

Since  we  have  not  presented  any  results  using  the  Newton-Raphson 
method  to  determine  the  Lagrange  multipliers,  we  will  not  discuss  the 
advantages  or  disadvantages  of  this  approach.  However,  in  using  this 
approach  the  coefficients  need  to  be  evaluated  for  all  the  elements 
and  constraints,  and  the  initial  values  of  the  Lagrange  multipliers  need 
be  assumed.  Also  in  the  Newton-Raphson  approach  it  is  generally  diffi¬ 
cult  to  eliminate  the  equations  associated  with  the  passive  constraints. 

The  last  approach,  based  on  using  simple  relations  to  determine  the 
Lagrange  multipliers,  is  computationally  efficient.  It  requires  less  CP 
time  to  perform  the  optimization  phase  of  the  algorithm.  With  this 
approach  it  is  possible  to  combine  all  the  virtual  loads  associated  with 
all  the  constraints  into  one  equivalent  virtual  load,  since  the  Lagrange 
multipliers  can  be  estimated  before  evaluating  the  coefficients  Q^-  or 
R.jj.  Particularly  for  the  stress-constrained  problem,  if  the  initial 
values  of  the  Lagrange  multipliers  are  assumed  to  be  proportional  to  the 
forces  in  the  bars  and  the  exponential  recurrence  relation  is  used  to 
modify  the  design  variables,  use  of  Eq.  132  leads  to  the  near-optimum- 
weight  design  with  the  least  computational  effort.  However,  if  the  con¬ 
straints  imposed  on  the  structure  are  mixed,  such  as  displacements  and 
stresses,  this  algorithm  needs  a  large  number  of  iterations  and  the  scaled 
weight  of  the  structure  may  violently  oscillate.  The  approximate  relation 
(Eq.  140)  to  determine  the  Lagrange  multiplier  for  the  FSD  design  gives 
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a  near-mi nimum-weight  design  if  the  minimum-weight  design  is  a  FSD.  The 
FSD  algorithm  (Eq.  116)  does  not  lead  to  the  correct  minimum-weight  design 
if  there  are  large  differences  in  the  maximum  allowable  stresses  in  the 
elements. 

In  the  optimization  of  a  structure  one  may  desire  to  obtain  a  design 
that  satisfies  the  theoretical  optimality  criterion.  However,  in  a  practical 
design  problem,  the  true  optimum  may  not  always  be  feasible  owing  to 
manufacturing  considerations.  In  many  practical  cases  it  may  be  useful 
to  obtain  a  design  that  is  a  near-minimum-weight  design  with  a  good  distri- 
bution  of  material  rather  than  a  design  satisfying  the  theoretical  optimality 
criterion.  For  many  problems  such  a  design  can  be  obtained  by  including 
only  a  small  number  of  constraints  into  the  potentially  active  category 
and  treating  the  other  constraints  as  passive  constraints.  A  good  example 
of  this  is  the  displacement-stress-constraint  problem.  For  such  a  problem 
it  is  advantageous  to  treat  all  the  stress  constraints  as  passive  con¬ 
straints  instead  of  including  them  in  the  constraint  equations  and  evaluat¬ 
ing  their  gradients  and  associated  Lagrange  multipliers.  Since  the  number 
of  active  stress  constraints  is  large  compared  to  the  displacement  con¬ 
straints,  this  simplification  saves  substantial  CP  time.  This  and  other 
simple  relations  proposed  in  the  previous  sections  to  determine  the 
Lagrange  multiplier  minimizes  the  total  computational  effort  needed  to 
design  an  optimum  structure. 
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APPENDIX 

THREE-BAR  TRUSS  PROBLEM 


A  three-bar  truss  shown  in  Figure  A-l  with  the  variation  of  loads 
and  constraints  will  be  used  to  illustrate  the  application  of  the 
optimality  criterion  approach  to  design  a  minimum-weight  structure.  The 
problem  is  small  and  most  of  the  required  calculations  can  be  done  by 
using  a  calculator. 


Figure  A-l.  Three-Bar  Truss 


Take  horizontal  and  vertical  displacements  of  the  node  P  are  denoted 

by  U  and  U  .  The  forces  applied  at  the  node  in  the  X  and  Y  directions 

x  y 

are  Fx  and  F  ,  respectively.  The  displacements  U1 ,  and  U2  and  U3  of  the 
point  P  in  the  directions  parallel  to  the  length  of  the  three  bars  are 
given  by 

U,  =  U  cos<p  +  U  cose 

u2  ■  u*  (A-l 

U.  .  -u?  cos$  +  U  cose 
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The  bar  forces  are  given  by 

A3EU3 

It,  2  l2  3  t3 

where  A^,  A2,  and  A3  are  the  cross-sectional  areas  of  the  three 
elements,  and  *2,  and  £3  are  their  respective  lengths. 


(A-2) 


The  stresses  in  the  elements  can  be  written  as 


The  equilibrium  equations  can  be  written  by  taking  a  summation  of 
the  forces  at  node  P.  This  gives 


T,  co s<|>  -  T_  cos<!>  =  F 

1  J  X  (A-4) 

T^  cose  +  T2  +  T3  cose  = 

Using  Equations  A-l  and  A-2,  these  equations  can  be  written  as 


h 

I 

IM 

I - 

ro 

1 

iv 

iFyl 

where 


2 

E  cos  <)> 


E  cos<t>  cose 


2  ^2 
E  cos  6  +  — 

*2 


E 


(A-6) 


Equations  A-5  are  the  load-displacement  relations  and  are  equivalent 
to  Equation  1  in  Section  I.  Solving  Equation  A-5,  the  nodal  displacements 
are  given  by 
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The  displacements  Ux 


and  U 

y 


are  also  be  given  by 


U..  = 


3 

E 

i=l 


Ti  s 

AiE 


u  = 


3 

E 

i=l 


Ti  *7 

AiE 


(A-8) 


1  2 

where  t.  and  t..  are  the  forces  in  the  ith  bar,  due  to  a  virtual  load  of 
unit  magnitude,  applied  at  node  P  in  the  x  and  y  directions,  respectively. 


The  equilibrium  equations  for  the  virtual  load  system  {f  ,  f  }  can  be 

x  y 

written  by  replacing  the  vector  {F  ,  F  }  in  Equation  A-5  by  {f  ,  f  }.  If 

*  y  x  y 

the  displacements  and  forces  due  to  the  virtual -load  system  are  denoted  by 
lower-case  letters,  then  the  expressions  for  these  quantities,  can  be 
written  by  replacing  the  capital  letters  by  the  lower-case  letters  for  the 
forces  and  displacements,  in  Equations  A-l ,  A-2  and  A-7. 

The  weight  of  the  three  bar  truss  is  given  by 


W  =  pJAjftj  +  A2*2  +  A3*3)  (A-9) 

The  solutions  of  the  problems  and  the  associated  algorithms  are 
discussed  in  subsequent  pages.  The  details  of  the  iterative  steps  are 
self  explanatory  and  follow  the  steps  given  in  Section  VIII. 


Problem  1.  A  Single-Displacement  Constraint 


E  =  10'  psi 


C]  =  0.05" 


I 


10,000  lb 

p  =  0.1  lb/in3  Am.n  =  0.1  in2 


°1  =  CT2  =  °3  =  ^5000  psi 
C1  =  Ux  or  Uy 


C,  <  C, 


•TV 
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This  problem  is  solved  by  using  four  different  algorithms.  All 
algorithms  lead  to  the  same  minimum-weight  design,  however,  the  number  of 
iterations  required  for  each  algorithm  is  not  the  same.  This  is  a 
single  constrained  problem  with  only  one  constraint  active  for  all 
iterations. 

Problem  1 (a). 

The  design  variables  are  modified  by  using  the  recurrence  relation 
(Equation  43). 


with  the  step-size  parameter  r  =  2.  The  Lagrange  multiplier  is  determined 
from  the  relation  (Equation  55), 


y  •'Qilpit1 

i=l 

^  =  —  TT 


Ci  -  C, 


i 

The  areas  A1 ,  A2  and  A3  at  the  beginning  of  each  iteration  are  given 
in  rows  1  through  3,  and  the  scaled  areas  are  given  in  rows  44  through  46. 
The  scaled  weight  after  each  iteration  is  given  in  row  48.  The  Lagrange 
|  multiplier  A 1  is  calculated  in  row  71. 


Row  48  shows  that  the  minimum  weight  is  113.7  lb  and  the  associated 

2  2 

cross-sectional  areas  are  A^  =  15.07  in  .  A^  =  0.1  in  .  and  A^  =  0.9318 

in^.  The  deflection  in  the  x  direction  at  node  P  is  the  active  constraint, 

and  A.  at  the  optimum  is  2141.0.  This  solution  satisfies  Equation  56 
1  *  _  * 
which  gives  the  relation  between  A^ ,  W  ^  ,  W  ,  and  C-j . 


I 


I 


I 


Problem  1(a)  Cont'd. 

CYCLE  t  CYCLE  t  CYCLE  S  CYCLE  »  CYCLE 


Problem  1(a)  Cont'd. 

CYCLE  1  CYCLE  S  CYCLE  3  CYCLE  %  CYCLE  *  CYCLE  6  CYCLE  T  CYCLE  • 


7$0tE»01t  .7504E»I1!  .7SB4EMQ1  .75*4E»0«t  .79B4E*«II 
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Problem  1(b). 

The  design  variables  lire  modified  by  using  the  recurrence  relation 
(Equation  43) 


1/r 

k 


with  the  step-size  parameter  r  =  2.  The  Lagrange  multiplier  is  determined 
from  the  relation  (Equation  57), 


The  arrangement  of  the  rows  is  the  same  as  for  Problem  1(a).  The  number 
of  iterations  required  to  obtain  the  minimum  weight  is  nearly  equal  to 
that  for  Problem  1(a). 
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1(b)  Cont'd. 


Problem  1(b)  Cont'd. 

cycle  1  cycle  *  cycle  j  cycle  *  cycle  *  cycle  6  ctcle  t  cycle 


CYCLE  i  CYCLE  *  CYCLE  3  CYCLE  8  CYCLE  %  CYCLE  «  CYCLE  7  CYCLE  8 
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Problem  1 (c) . 

The  design  variables  are  modified  by  using  the  recurrence  relation 
(Equation  45), 


a?+’  =  a(  ('  *Hx.  .? 


Vi‘i 


with  the  step  size  parameter  r  =  2.  The  Lagrange  multiplier  is  determined 
from  the  relation  (Equation  57) 


x,  = 


C1  '  C1 


1  ?  Qfi 


1=1  VA 


The  convergence  is  slower  than  for  Problems  1(a)  or  1(b).  After  eight 
iterations  the  minimum  weight  is  113.7  lb,  but  comparing  the  areas  of  the 
members  with  Problem  1(a)  shows  that  the  solution  has  not  as  yet  converged. 
An  additional  2  or  3  iterations  are  needed. 
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Problem  1(c)  Cont'd. 

CYCLE  1  CYCLE  t  CYCLE  3  CYCLE  t 
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The  design  variables  are  modified  by  using  the  recurrence  relation 
(Equation  158), 


with  the  step-size  parameter  r  =  2.  The  Lagrange  multiplier  is  determined 
from  the  relation  (Equation  57) 


X 


1 


* 


1=1  P.4. A? 

l  l  l 


This  algorithm  is  equivalent  to  the  algorithm  based  on  the  use  of  the 
reciprocal  design  variable.  The  convergence  for  this  algorithm  is  slower 
than  for  all  the  three  previous  problems.  After  eight  iterations  the 
minimum  weight  is  123.0  lb  and  needs  additional  iterations  for  the  solution 
to  coverge  to  the  minimum. 


AFWAL-TR -81-31 24 


CYCLE  «  CYCLE  6  CYCLE  T  CYCLE  • 
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Problem  2.  Multiple  Displacement  Constraints. 


20,000  lb 


E  =  107  psi  p  =  0.1  lb/in3  Amin  =  0.1  in2 


°1  =  °2  =  03  =  25,000  psi. 

=  0.05  in.  C  2  =  0.05  in. 


C1  -  C1 


This  problem  is  solved  by  using  three  different  algorithms.  All 
algorithms  lead  nearly  to  the  same  minimum-weight  design,  but  the 
convergence  behavior  is  not  the  same.  The  two  constraints  are  active 
throughout  the  iterative  history. 
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S 


i 


Problem  2(a). 

The  design  variables  are  modified  by  using  the  recurrence  relation 
(Equation  71 ) 


A1?*1  =  Aj  (  1  +  -  (  Z  A .  -1 

'  ,l  rVj=l  VlN  '/k 


with  the  step-size  parameter  r  =  2.  The  Lagrange  multipliers  are  determined 
from  the  relation  (Equation  86) 


— 1 
> 

B 

i 

Ri 

1 — 
CO 

o 

L 

l  x2  ) 

•  R2. 

where 


A  =  E 


Ojl  ‘ill 
p.k.A,3 


(Active  Elements) 


B  =  E 


Qil  qi2 


(Active  Elements) 


C  =  E 


Qi2  Qi 2 
pi£iAi 


(Active  Elements) 


R-j  =  2(C]  +  B1 


C1 )  +  B] 


Ro  = 


2(C2  +  B2 


c2)  +  b2 


★  w-i 

C1  =  E  A“ 


(Passive) 
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*  ^i2 

C2  =  £  A“ 

(Passive) 

B  -  s  !!ll 

B1  1  A. 

(Active) 

^i2 

b2  *  £  a“ 

(Active) 

The  Lagrange  multipliers  are  given  by 


x 


1 


B-R2  -  C-R, 
B2  -  AB 


where 


B*  R-j  -  A*R2 

- 2 - 

B  -  AB 


Ml 


j 

_  ill 


and  Q 


ViS 


12 


1  2 

and  ti  are  the  forces  in  the  bars  dur  to  a  unit  load  applied  at  node  P 
in  the  X  and  Y  direction  respectively. 

The  weight  of  the  structure  after  eight  iterations  was  151.1  lb. 

One  more  iteration  was  needed  to  reach  the  optimum  design.  The  optimum 
weight  of  the  structure  was  150.7  lb  with  A1  =  14.14  in2,^,A2  =  10.10  in2 
and  A^  =  0.1  in2.  For  the  optimum  design,  x^  =  1000  and  x2  =  2000.  This 
design  satisfies  the  relation 


W-W* 


m 


E 

j=l 


which  is  valid  at  optimum  solution. 


Problem  2(a)  Cont‘d. 

crcte  t  crcte  *  crcte  I  Crete  t  crcte  *  Crete  6  Crete  r  Crete  • 
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Problem  2(a)  Cont'd. 

CYCLE  X  CYCLE  t  CYCLE  S  CYCLE  t  CYCLE  9  CYCLE  6  CYCLE  ?  CYCLE  I 
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AFWAL-TR-81-3124 


125 


ltME»01t  .122OE»01t  .130&£*01t  .1357E*0H  .13R«E«01( 


AFWAL -TR -81-31 24 


Problem  2(b). 


The  design  variables  are  modified  by  using  the  relation  (Equation  158). 


£ 

l 


i  /  m  Q. . 

'  -  \(  *  >A 

,  \J=1 


with  the  step-size  parameter  r  =  2.  The  Lagrange  multipliers  are 
determined  by  using  the  same  equations  as  for  problem  2(a)  (See 
Equations  84  and  161).  From  row  48  it  is  seen  that  the  weight  of  the 
structure  after  eight  iterations  is  156.0  lb,  and  a  few  more  iterations 
are  required  in  order  to  reach  the  minimum  weight  of  150.7  lb.  The 
algorithm  used  for  this  problem  is  equivalent  to  the  one  derived  by  using 
the  reciprocal  design  variable.  This  algorithm  is  generally  slower  to 
converge  than  the  one  used  for  problem  2(a). 
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Problem  2(b)  Cont'd. 

CYCLE  1  CYCLE  t  CYCLE  3  CYCLE  t  CYCLE  9  CYCLE  6  CYCLE  T  CYCLE  • 


Problem  2(b)  Cont'd. 
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Problem  2(c). 

In  this  problem  the  relative  design  vector  cu  (See  Equation  17)  is 
used.  This  is  modified  by  using  the  relation 


k+1  k  /  Qil  \ 

°1  "  °i  l  2  ) 


1/2 


i-i  1  '  k 

and  the  Lagrange  multipliers  are  estimated  from  the  relation  (Equation  74) 


For  the  first  iteration  X-|  =  X2  =  1.0.  The  Lagrange  multipliers  and  the 
relative  design  variables  o^.  are  normalized  after  each  modification  so 
that  the  maximum  a.  or  X.  is  equal  to  unity.  A  single  virtual -load  vector 
with  f  =  x. .  1  and  f  =  x?.l  is  used. 

A  I  j  t 

After  eight  iterations  the  weight  of  the  structure  is  150.8  lb 

? 

(See  Row  No.  48),  but  A,  for  this  design  is  0.0999  in  ,  which  should  be 
2  J 

equal  to  0.1  in  .  The  minimum-weight  design  is  nearly  the  same  as  the 
one  for  Problem  2(a). 


Problem  2(c)  Cont’d. 

CYCLE  1  CYCLE  t  CYCLE  3  CYCLE  fc  CYCLE  3  CYCLE  «  CYCLE  T  CYCLE  • 
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Problem  3.  Stress  Constrained  Problem. 


o i  =  10,000  psi  =  40,000  psi 


a ^  =  10,000  psi 


This  problem  is  solved  by  using  three  different  algorithms.  First 
the  structure  is  designed  by  using  the  F.S.D.  method,  then  the  other  two 
approaches  are  used.  For  this  problem  the  F.S.D.  algorithm  does  not  lead 
to  the  real  minimum-weight  design. 
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Problem  3(a). 

The  design  variables  are  modified  by  (Equation  117) 


This  is  the  Fully  Stressed  Design  algorithm.  The  weight  of  the 
structure  with  the  same  size  for  all  the  elements  is  247.5  lb.  After  the 
first  iteration  the  weight  reduced  to  158.4  lb.  In  the  second  iteration 
the  weight  increased  to  246.7  lb  and  then  decreased  with  subsequent 
iterations  until  it  reached  198.5  lb.  The  cross-sectional  areas  for  the 
design  with  a  weight  equal  to  198.5  lb  are  =  21.07  in2,  A?  =  0.1  in2 
and  A,  =  6.92  in2.  The  stress  in  the  second  element  is  equal  to  20,000  ps 

J  /i 

in  the  other  two  elements,  10,000  psi .  The  size  of  A^  is  dictated  by  the 
minimum-size  requirement.  This  is  not  an  optimum  design  (See  Problem 
3(b)  or  3(c)). 
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Problem  3(a)  Cont‘d. 
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Problem  3(b). 

In  this  problem  the  relative  design  variable  is  used.  It  is 
modified  by  using  the  relation 


k+1 

flj 


1/2 

k 


The  Lagrange  multipliers  are  estimated  by  using  the  relation  (Equation  132) 


3=1 ,2,3 


For  the  first  iteration  the  A.'s  are  assumed  to  be  proportional  to  the 

J 

forces  in  the  bars.  The  relative  design  variables  and  the  Ligrange 
multipliers  are  normalized  after  each  iteration  so  that  the  maximum 
or  A.  is  equal  to  unity.  A  single  virtual  load  vector  is  used  with 


AFWAL-TR-81  -31 24 


T. 

The  ratio  is  equal  to  +1  or  -1  depending  on  whether  the  force  in 

the  bar  is  Tension  or  Compression. 

Row  48  shows  the  minimum  weight  of  the  structure  to  be  128.6  lb.  This 
is  achieved  with  five  iterations.  With  additional  iterations  the  weight 
increases.  The  cross-sectional  areas  for  the  minimum-weight  design  are 

2  2  2 

A^  =  14.22  in  ;  A2  =  5.468  in  ;  A^  =  0.1003  in  .  The  Lagrange  multipliers 
associated  with  this  design  are  =  1.0,  =  0.086  and  a3  =  0.01.  For 

convenience  the  flexibility  coefficient  is  denoted  by  Qij  instead  of  Rij. 


Problem  3(b)Jtont'd. 

CYCLE  1  CYCLE  E  CYCLE  3  CYCLE  %  CYCLE  9  CYCLE  6  CYCLE  7  CYCLE  • 

IV  HE  *111  .10MEH1I  .lMOEMlt  .10V0E»01t  .1I00E»0U  .1MOEH1I  .10HE*»1I  .1HBE«V1I 
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Problem  3(b)  Cont'd 

CYCLE  1  CYCLE  *  CYCLE  s  CYCLE  % 
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M! 


Probjem  3(b)  Cont'd. 

CYCLE  ft  CYCLE  t  CYCLE  3  CYCLE  t  CYCLE  9  CYCLE  •  CYCLE  7  CYCLE  • 

r«rftc*Mt  •••/•E-flu  .ms’iii  .tttrMii  .inu«in  .iehmh  .jeiye-mi  .mkhii 
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Problem  3(b)  Cont'd. 

me  1  CYCLE  t  CYCLE  3  CYCLE  ■»  CYCLE  3  CYCLE  6  CYCLE  7 


Problem  3(b)  Cont'd 
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Problem  3(c) . 


In  this  problem  the  design  variables  are  modified  by  using  the 
relation  (Equation  131) 


A 


-1 


k 


with  step  size  parameter  r  =  2.  The  equations  used  for  evaluating  the 
Lagrange  multipliers  are  the  same  as  those  of  Problem  2(a).  The  stress 
constraints  in  elements  1  and  3  are  assumed  to  be  potentially  active  in 
order  to  reduce  the  size  of  the  problem  solution.  The  stress  in  element 
2  is  inactive,  and  assuming  it  to  be  potentially  active  would  only 
increase  the  size  of  the  problem.  For  element  1  the  virtual  load  is 


and  for  element  3  the  virtual  load  is 


The  flexibility  coefficients  are  denoted  by  Qij  instead  of  Rij  for 
convenience.  The  weight  of  the  optimum  design  is  126.06  lb,  and  the 
cross-sectional  areas  of  the  three  members  are  A<  =  14.24  in  ,  A?  *  4.929 

0  0  x  c 

1nc  and  Ag  =  0.1  in  .  The  stresses  in  the  three  elements  are  10,000 
psl,  20,000  psi  and  10,000  psi  respectively.  The  stress  constraint  In 
element  1  is  active  for  all  iterations,  but  the  stress  constraint  in 
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element  3  becomes  active  only  after  the  fourth  iteration.  The  stress  in 
element  2  is  not  active. 

Comparing  this  design  with  the  one  obtained  by  using  the  FSD 
algorithm,  one  finds  that  the  two  designs  are  not  the  same  although  the 
stress  distribution  for  them  is  the  same.  The  design  obtained  by  the  FSD 
algorithm  is  non-optimum. 


a 
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